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Abstract. We prove duality relations for two interacting particle systems: the g-deformed to- 
tally asymmetric simple exclusion process (q-TASEP) and the asymmetric simple exclusion process 
(ASEP). Expectations of the duality functionals correspond to certain joint moments of particle 
locations or integrated currents, respectively. Duality implies that they solve systems of ODEs. 
These systems are integrable and for particular step and half stationary initial data we use a nested 
contour integral ansatz to provide explicit formulas for the systems' solutions and hence also the 
moments. 

We form Laplace transform like generating functions of these moments and via residue calculus 
we compute two different types of Fredholm determinant formulas for such generating functions. For 
ASEP, the first type of formula is new and readily lends itself to asymptotic analysis (as necessary 
to reprove GUE Tracy- Widom distribution fluctuations for ASEP), while the second type of formula 
is recognizable as closely related to Tracy and Widom's ASEP formula [3n]-[12]. For q-TASEP both 
formulas coincide with those computed via Borodin and Corwin's Macdonald processes [8]. 

Both g-TASEP and ASEP have limit transitions to the free energy of the continuum directed 
polymer, the logarithm of the solution of the stochastic heat equation, or the Hopf-Cole solution 
to the Kardar-Parisi-Zhang equation. Thus, g-TASEP and ASEP are integrable discretizations of 
these continuum objects; the systems of ODEs associated to their dualities are deformed discrete 
quantum delta Bose gases; and the procedure through which we pass from expectations of their 
duality functionals to characterizing generating functions is a rigorous version of the replica trick in 
physics. 
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1. Introduction 

One dimensional driven diffusive systems play an important role in both physics and mathematics 
(see e.g. [Ml IMl E])- As physical models they are used to study mass transport, traffic flow, 
queueing behavior, driven lattice gases, and turbulence. Their integrated current defines height 
functions which model one-dimensional interface growth. In certain cases they can be mapped into 
models for directed polymers in random media and propagation of mass in a disordered environment. 



1 



FROM DUALITY TO DETERMINANTS FOR Q-TASEP AND ASEP 



2 



The particle systems provide efficient means to implement simulations of these various types of 
systems and, in some rare cases, yield themselves to exact and rigorous mathematical analysis. 

This article is concerned with two interacting particle systems - g-TASEP with general particle 
jump rate parameters, and ASEP with general bond jump rate parameters - which contain rich 
mathematical structure. Presently we seek to shed light on structure which exists in parallel for both 
of these systems. We demonstrate duality relations (see Definition I2.ip for both of these systems 
directly from their Markovian dynamics: g-TASEP is dual to a totally asymmetric zero range 
process TAZRP (Theorem 12. 2p whereas ASEP is self dual (Theorems 14.11 and 14. 2p . A consequence 
of duality is that expectations of a large class of natural observables of these systems evolve according 
to systems of ODEs. 

For g-TASEP, the duality result is, to our knowledge, new. When all particle jump rate parameters 
are equal, dynamics of g-TASEP can be encoded via a quantum integrable system in terms of q- 
Bosons [Bl [34]. For ASEP with all bond jump rate parameters equal, the ASEP self duality was 
observed by Schiitz [35] and further developed in [19] via a spin chain representation of ASEP (which 
is related to the XXZ model - a well studied quantum integrable system). Our results apply for 
general rates and proceed directly via the Markov dynamics. 

The most surprising observation of this article is that, for certain initial data called step and half 
stationary, we are able to explicitly solve the systems of ODEs for g-TASEP and ASEP in terms 
of simple nested-contour integrals. For g-TASEP this works for the full generality of particle jump 
rate parameters, whereas for ASEP we must assume all bond jump rate parameters to be equal 
at this stage and henceforth. For g-TASEP, the integral representations of the solution can also 
be obtained via the formalism of Macdonald processes [HJ [10] , while for ASEP we were guided by 
analogy and results of [T9] . 

Let us state the simplest versions of these formulas, focusing just on step initial data. We also 
informally introduce the dynamics of g-TASEP and ASEP. 

The (7-TASEP is a continuous time Markov process x{t). Particles occupy sites of Z and the 
location of particle i at time t is written as Xi{t) and particles are ordered so that Xi{t) > Xj{t) for 
i < j. The rate at which the value of Xj(t) increase by one (i.e. the particle jumps right by one) is 
aj(l — qi^i-i(*)-^i(*)-i); all jumps occur independently of each other according to exponential clocks. 
Here q E [0, 1) represents the strength of the repulsion particle Xi feels from particle Xi-i. For the 
purpose of this introduction we restrict to = 1 and consider only step initial data where particles 
start at every negative integer location and nowhere else (i.e., for i > 1, Xi{0) = —i). The following 
result appears as Corollary 12.121 

Theorem 1.1. Consider q-TASEP with step initial data and particle jump rate parameters ai = 1. 
Then for any k > 1 and ^^l > ^2 > • • • > > 0, 



E 



n 



(-l)*^g'*V f f ^ ZA-ZB-^,^ ^ ^_n^^{q-l)tz,dZj 



/■■■/ n If^na-,)- 

l<A<B<k ^ ^ ^ j=l 

where the integration contour for za contains {qzB}B>A md 1 hut not 0. 

The ASEP (occupation process) is a continuous time Markov process r]{t) = {rjx{t)}x&i- The 
r]x (t) are called occupation variables and are 1 or based on whether there is a particle or hole at x 
at time t. The dynamics of this process is specified by non-negative real numbers p < q (normalized 
by P + Q = 1) and uniformly bounded (from infinity and zero) rate parameters {a^jxez- For each 
pair of neighboring sites {y, y + 1) the following exchanges happen in continuous time: 

r] i-> r/^'S^+i at rate Oyp if (??y, r/^+i) = (1,0), 

r] h-^ 7y!^'f+i at rate ayq if {r]y,r]y+i) = (0, 1), 
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where r/^'2'+^ denotes the state in which the value of the occupation variables at site y and y + 1 
are switched, and all other variables remain unchanged. All exchanges occur independently of each 
other according to exponential clocks. For the purpose of this introduction we restrict to a^: = 1 
and consider only step initial datsQ where r}x{Q) = lx>i- Assume < p < q and let t = p/ q < \. 
Finally, let N^it) = ^y^^ Vyii) record the number of particles to the left of position x + 1 at time t. 
The following result on ASEP appears as Theorem 14.191 

Theorem 1.2. Consider ASEP with step initial data and all bond rate parameters = 1. Then 
for all n> 1 and x G Z, 



E 



n{n~l) 
T 2 



1 



n 



ZA - ZB 



n 



l<A<B<7i ^ " i=l 



exp 



Zi{p - qf 



{zi + l){p + qzi) 



l + Zj 
1 + Zi/T 



dzi 



where the integration contour for za includes 0, — r hut does not include —1, or {tzb}b>a (see 
Figure\^for an illustration of such contours). 

These expectations contain sufficient information to uniquely characterize the distribution of the 
location of a given collection of particles (after the system has evolved for some time) in each of 
these systems. Focusing on a single Xn{t) or Nx{t) distribution we can concisely characterize this 
via generation functions of suitable expectations. There are two types of generating functions we 
consider - both related to g-deformed (or for ASEP r-deformed) Laplace transform introduced by 
Hahn [I?] in 1949. 

These generating functions are naturally suggested from the nested structure of the contour 
integral formulas for these expectations. There are two ways to deform the nested contour integrals 
so all contours coincide. Accounting for the residues encountered during these deformations, we are 
led to two types of formulas for expectations: Those involving partition-indexed sums of contour 
integrals and those involving sums of contour integrals indexed by natural numbers. 

Using the partition-indexed formulas we prove that the first generating function is equal to a 
Fredholm determinant which we call Mellin-Barnes type. The following result is contained in 
Theorem [ 



Theorem 1.3. Consider ASEP with step initial data and all bond rate parameters = 1. Then 
for all X eZ and C &C\R+ 



E 



1 



(CT^-W;r) 



det(/ + 



ASEP\ 



where (a;r)oo = (1 — a)(l — ra) • • • , and where the space on which K^^^^ acts can be found in 
the statement of Theorem \5.3[ The operator K^^ is defined in terms of its integral kernel 



Kf^P{w,w') = ^ [ T{-s)r{l + s){-Cy 



9x{w) 



-ds. 



gxiT^w) w' — T^W 

(The contour Dji d is specified in the statement of Theorem \5.S\) . The function gx{z) is given by 



gx[z) = exp 



{q-p)t 



Z + T 



Z + T 



Observe that the step initial data for g-TASEP involves particles to the left of the origin, whereas for ASEP it involves 
particles to the right of the origin. We decided to keep these conventions to be consistent with previous works on the 
subject. 
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This type of formula lends itself to rigorous asymptotic analysis. For ASEP this formula is new 
and in Section [9] we demonstrate how it can be used to recover Tracy and Widom's celebrated 
fluctuation result [31] which states that 



Here j = q — p is assumed to be strictly positive and -Fgue is the GUE Tracy- Widom distribution. 
The case when 7 = 1 (g = 1 and p = 0) was proved earlier by Johansson |21j . Theorem 11.31 also 
allows to access (under a certain weakly asymmetric scaling) the narrow wedge KPZ equation one 
point formula [U [33] . 

For g-TASEP such a Mellin-Barnes type formula was obtained from the theory of Macdonald 
processes [8], Theorem 4.1.40. It should be possible to use this Fredholm determinant to prove 
cube-root GUE Tracy- Widom fluctuations for the current past the origin in g-TASEP. This has not 
yet been done, though in an stationary version of the TAZRP associated to g-TASEP gaps, the cube- 
root fluctuation scale is shown in [3] (via a different approach) . In [HJ [9] the q-TASEP Mellin-Barnes 
type Fredholm determinant formula is used (via a limit transition) to write the Laplace transform 
of the O'Connell-Yor semi discrete polymer partition function [29]. Then [3 |9] perform rigorous 
asymptotic analysis to show cube-root GUE Tracy- Widom free energy fluctuations as well as to 
provide a second rigorous derivation of the narrow wedge KPZ equation formula (first rigorously 
derived in [T]). From the perspective of asymptotics, this second approach is a little less involved 
than that of [T]. 

On the other hand, using the formulas of the second type (i.e. deforming contours in Theorems ll.il 
and [L2] differently) we prove that the second generating function is equal to a Fredholm determinant 
which we call Cauchy type. The ASEP Fredholm determinant Tracy and Widom derived in [39^ BO] 
is also of this type (and in fact, after inverting the Cr-Laplace transform we recover the same 
formula as in [39^ I40j). Asymptotic analysis of this type of determinant is not as straightforward as 
the Mellin-Barnes type. In [H], Tracy and Widom employ a significant amount of post-processing 
to turn this type of formula into one for which they could perform asymptotic analysis. The final 
formula still involves a complicated term related to the Ramanujan summation formula (as observed 
in [33]). One should note that while we do recover (among other formulas) the Tracy- Widom ASEP 
Fredholm determinant formula, our approach via duality is entirely different, our contour integral 
ansatz is not a version of the coordinate Bethe ansatz and, along the way, we gain access to other 
information about ASEP, like joint moment formulas. The Cauchy type Fredholm determinant 
formula for g-TASEP was also first derived in [8] via Macdonald processes. 

In short, by utilizing duality for g-TASEP and ASEP, we are able to provide a short and di- 
rect route from Markov dynamics to Fredholm determinant formulas characterizing single particle 
location or single integrated current distributions. 

Both g-TASEP and ASEP are integrable discretizations of the KPZ equation. As stochastic 
processes they converge to the Hopf-Cole solution to the KPZ equation [5l[Il[27]. The systems of 
ODEs associated with their duality have limit transitions to the attractive quantum delta Bose gas 
which describes the evolution of joint moments of the stochastic heat equation (whose logarithm is 
the KPZ equation and which describes the partition function for the continuum random polymer). 

An advanced version of the popular physics polymer replica trick attempts to recover the Laplace 
transform of the one point distribution of the solution to the stochastic heat equation in terms 
of its moments (see Section [6.4p . However, the moments grow far too quickly to characterize this 
distribution and hence drawing conclusions from them is mathematically unjustifiable and in any 
case, risky. Nevertheless, Dotsenko [15] and Calabrese, Le Doussal and Rosso [TT] were eventually 




(1) 



FROM DUALITY TO DETERMINANTS FOR Q-TASEP AND ASEP 



5 



able to use this trick to recover the exact formulas of [Tl[33]. Given the results of our work, the non- 
rigorous replica trick manipulations can be seen as shadows of the rigorous duality to determinant 
approach developed presently for both g-TASEP and ASEP. That is to say that by going to a 
suitable discrete approximation we are able to rigorously recover analogs of Laplace transforms 
from moments and then in the limit transition these converge to formulas for the stochastic heat 
equation's Laplace transform. The replica trick has proved computationally useful (see for instance 
[20j). thus providing additional motivation for the present work. 

The limit transition of g-TASEP to the O'Connell-Yor semi-discrete directed polymer [29] (and 
associated semi-discrete stochastic heat equation) is explored in Section [6l Under that limit transi- 
tion, duality becomes the replica approach and the duality system of ODEs become a semi-discrete 
version of the delta Bose gas. The nested contour integral ansatz provides means to succinctly 
compute the solution to the Bose gas. The Fredholm determinants for q-TASEP limit to Fredholm 
determinants for the Laplace transform of the polymer partition function. 

While a variety of probabilistic systems arise as degenerations of Macdonald processes, ASEP is 
not known to be one of them. For ASEP, it is not known what, if anything, replaces this additional 
integrable structure endowed to g-TASEP from its connection to symmetric functions. However, it 
is compelling that both g-TASEP and ASEP have duality relations and that the associated systems 
of ODEs can both be solved via a nested contour integral ansatz. This leads one to ask whether 
g-TASEP and ASEP can be unified via a theory even higher than Macdonald processes. Spohn [33 
has coined the term stochastic integrability to describe stochastic processes which display a great 
deal of integrable structure. Both g-TASEP and ASEP are clear examples of such processes and 
the contributions of this work provide an additional layer to that integrability. 

1.1. Outline. The paper is organized as follows. In Section [2] we prove duality for g-TASEP and 
explicitly solve the associated systems of ODEs via a nested contour integral ansatz. In Section[3]we 
provide a general scheme to go from such nested contour integral formulas to two types of Fredholm 
determinants and apply this to Theorem 11.11 In Section [4] we prove duality and nested contour 
integral formulas for ASEP. In Section [5] we explain the passage from Theorem 11.21 to Theorem 
11.31 Section [6] deals with a degeneration of g-TASEP to a semi-discrete directed polymer. Section 
[7] collects necessary combinatorial facts. Section [8] proves a uniqueness result for the system of 
ODEs associated with ASEP duality. Section [9] provides critical point analysis of the Fredholm 
determinant in Theorem 11.31 as necessary to obtain ([1]). 

1.2. Notations. We fixed a few notations used throughout this paper. The imaginary unit l = 



The indicator function of an event E is denoted by either 6e or 1e- We write = 1 if 
ax = 1 for all x. All contours we consider are simple, smooth, closed and counter-clockwise oriented 
(unless otherwise specified). For a contour C we write aC as the dilation of C by a factor of a > 0. 
When we write that the integration contour for contains {qzB}B>A we mean that the contour 
contains the image of the zb contour dilated by q. Containment is strict so that if C contains a 
point a, then C separates a from infinity and the distance from C to a is strictly positive. 

1.3. Acknowledgements. We wish to thank P. Ferrari, T. Imamura, J. Quastel, C. Tracy and 
H. Widom for multiple discussions on the subject. Some work on this subject occurred during the 
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American Institute of Mathematics workshop on the Kardar-Parisi-Zhang Equation and Universality 
Class, and the Institute for Mathematical Sciences, NUS, workshop on Polymers and Related Topics. 
AB was partially supported by the NSF grant DMS-1056390. IC was partially supported by the NSF 
through PIRE grant OISE-07-30136 and DMS-1208998 as weh as by Microsoft Research through 
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2. Duality and the nested contour integral ansatz for g-TASEP 

The ^-deformed totally asymmetric simple exclusion process ((^-TASEP) is a continuous time, 
discrete space interacting particle system x{t). Particles occupy sites of Z and the location of 
particle i at time t is written as Xi{t) and particles are ordered so that Xi{t) > Xj{t) for i < j. 
The rate at which the value of Xi{t) increase by one (i.e. the particle jumps right by one) is 
aj(l — gi^i-i(*)-^i{*)-i); all jumps occur independently of each other according to exponential clocks. 
Here q € [0, 1), a,, > is particle i's jump rate parameter, Xi-i{t) — Xi[t) — 1 is the number of empty 
sites to its right (before particle Xi-i) and all jumps occur independently of each other. We will 
use and to denote expectation and probability (respectively) of the Markov dynamics with 
initial data x. When the initial data is itself random, we write E and P to denote expectation and 
probability (respectively) of the Markov dynamics as well as the initial data. We also use E and P 
when the initial data is otherwise specified. 

We presently focus on g-TASEP with particles xi > X2 > ■■■ > x^r. However, to ease the 
statement of results we include a virtual particle xo{t) = oo and define our state space as 

X'^ = {x = (xo, xi, . . . , xtv) G {oo} X :00 = Xo>Xi>--->XAr}. 

In this case, the dynamics are easily seen to be well-defined. Observe that the evolution of the 
right-most M < N particles performs q'-TASEP with M particles (i.e., particles are unaffected by 
those to their left). On account of this, it is easy to extend the dynamics to an infinite number of 
particles labeled xi > X2 > • • • (i.e., there is a right-most particle). When studying these infinite 
systems, it is generally enough to study related finite systems. 

For g-TASEP with N particles, the generator of x{t) acts on functions / : — t- M and is given 

by 

N 

(^,-tasep^)(^) = Y^ai{l - g--i-»-i) (/(f+) - /(f)) , (2) 

i=l 

where xf indicates to increase the value of Xj by one. Note that one may also write down a generator 
in terms of occupation variables and (as in [8]) show that for any initial data g-TASEP is, in fact, 
well-defined. 

The totally asymmetric zero range process (TAZRP) on an interval {0,1,..., A^} with site- 
dependent rate functions gi : Z>q — >• [0, oo) (with gi{0) = fixed) is a Markov process y{t) with 
state space 

y^ = (Z>o)^°'^'-'^>. 

The dynamics of TAZRP are given as follows: for each i G {1, • • • ,N}, yi{t) decreases by one and 
yi-i(t) increase by one (simultaneously) in continuous time at rate given by gi{yi(t)); for different 
i's these changes occur independently. Note that no particles leave site 0. The rate functions we 
consider are given by gi{k) = aj(l — q^). When all = 1, this model was first introduced in |34j . 
The generator of y{t) acts on functions h : Y'^ — >• R and is given by 

N 

(L«-TAZRP/,)(y-) = ^a,(l - qy^) {h{y^''-') - h{y)) (3) 

i=l 

where y*'*"^ indicates to decrease yi by one and increase yi-i by one. 

Observe that the gaps yi{t) = Xi(t) — Xi+i(t) — 1 of (7-TASEP evolve according to a TAZRP, but 
with boundary conditions that yo{t) = yN(t) = 00 for all t £ M+. Our work will not draw on this 
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xq X5 X4 X3 X2 xi Ug 2/5 2/4 ys 2/2 J/1 2/0 

Figure 1. Left: g-TASEP with six particles. The indicated jump of 2:4 occurs at 
rate 04(1 — g^) since the gap X3 — X4 — 1 = 2. Right: The dual TAZRP with sites 
{0, 1, . . . , 6}. The indicated jump occurs at rate 02(1 — q"^) since y2 = 2. 



obvious coupling. Rather, our statement of duality will provide a different relationship between x{t) 
and an independent y{t). 

2.1. Duality. Recall the general definition of duality given in Definition 3.1 of |24j . 

Definition 2.1. Suppose x{t) and y{t) are independent Markov processes with state spaces X and 
Y respectively, and let H{x, y) be a bounded measurable function on X xY . The processes x{t) 
and y{t) are said to be dual to one another with respect to H if 

¥.-[H{x{t),y)]=¥y[H{x,ym (4) 

for all X G X and y ^Y . Here W refers to the process x{t) started with x(0) = x (likewise for y). 

Theorem 2.2. The q-TASEP x{t) with state space X^ and particle jump rate parameters at > 0, 
and the TAZRP y{t) with state space Y^ and rate functions gi{k) = aj(l — g*^) are dual with respect 
to 

N 

Remark 2.3. The definition of H{x, y) means that = if yo > and if(x, y) = n,=i if 
2/0 = 0. 

Before proving Theorem 12.21 we define the following system of ODEs. 

Definition 2.4. We say that y) : M_|_ x Y^ — )■ M solves the true evolution equation with initial 
data ho{y) if 

(1) For all yeY^ and t G M+, 

(2) For ah yeY^ such that yo > 0, h{t; y) = for all t G M+; 

(3) For ally^Y^,h{0;y) = ho{y). 

Remark 2.5. The existence and uniqueness of global solutions to (A) is assured since it reduces 
to a finite system of linear ODEs, from which the result follows from standard methods |12j . 



Proof of Theorem \2.S\ . We claim first that for x and y fixed, 

^.-TASEP^(^^^ = L'^-TAZRP^(x, y), (5) 

where in the above expression, the generator on the left acts in the x variables and the generator 
on the right in the y variables. 
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To prove the claim is easy. Observe that 

N / N 

^<7-TASEP^(^^ y) = g^-i-^'-l) {qV' - 1) JJ 

i=l \ j=0 

N 

i=l 

Given the claim we may now check that [H{x{t),y)] and [H{x,y{t))] both satisfy the true 
evolution equation given in Definition 12. 4[ By the uniqueness of Remark l2.5l this implies the desired 
equality to complete our proof. That E^ [H{x,y{t))] satisfies this evolution equation follows from 
the definition of the generator of y{t). 

On the other hand, 

^E-^[H{xit),y)] = L'?-TASEP]E-[^(^(i)^y-)] 

= E^[L''-^^^^^H{x{t),y)] 
= E^"[L^-TAZRPi7(f(t),y)] 

= L'i'^^^^^E^[H{x{t),y)]. 

The equality of the first line is from the definition of the generator of x{t); the equality between the 
first and second lines is from the commutativity of the generator with the Markov semi-group; the 
equality between the second and third lines is from applying equality ([5]) to the expression inside 
the expectation; the final equality is from the fact that the generator ^^-TAZRP j-^^^ g^^^-g ^.j^g ^ 
coordinate and the expectation acts on the x coordinate. This shows that E^ [H{x{t), y)] solves the 
system of ODEs in the true evolution equation (checking the boundary condition and initial data 
is easy). □ 

2.2. Systems of ODEs. As a result of duality, we provide three different systems of ODEs to 
characterize [H{x{t),y)]. It is convenient to introduce an alternative way to write a TAZRP 
state y ^Y. 

Definition 2.6. For A; > 1, define 

W>o = {n= (ni,n2,...,nfc) S (l^o)'' : ni > n2 > . . . > > O}. 

For k = X^^il/j, associate to y E a vector n = n(y) G W^q where for i G {!,..., A^}, 
\{nj : Uj = i}\ = yi. Likewise, to a vector n G W^q we may associate y = y{n) G Y^ with the 
convention that yo{n) = 0. For instance, if A^ = 3, yi = 2, y2 = and 7/3 = 1 then /c = 3, ni = 3 and 
n2 = = 1. A vector n naturally splits into clusters, which are maximal groupings of consecutive 
equal valued elements. For instance if n = (4, 4, 2, 1) we would say there are three clusters with the 
cluster of 4 containing two elements, and the clusters of 2 and 1 containing only one elements each. 

Also, define the difference operator V/(n) = /(n — 1) — f{n). For a function f{n), Vj acts as V 
on the Hi variable. Finally let = (ni, . . . ,ni — 1, . . . , n^). 

Proposition 2.7. Let x G and x{t) be the q-TASEP started from f(0) = x. 

(A) True evolution equation: Ifh{t; y) : M_|_ x — > M solves the system of ODEs given in Definition 
\2.4\ with initial data ho{y) = H{x,y), then for all y G Y^ , E^ lH{x{t),y)] = h{t;y). 
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(B) Free evolution equation with boundary condition: If u : M_|- x (Z>o) 



solves: 



(1) For all n £ (Z>o)'' and t G 



d ^ 
— n) = (1 - 9) ^ an,Viu{t; n); 



(2) For all n G {Z 



>0) 



such that for some i € {1, . . . ,k — 1}, Ui 
Viu{t; n) = qVi+iu{t; n); 
(3) For all n G (Z>o)'^ such that = 0, u(t; n) = for all t G\ 



(4) For all n £ W!;q, n(0; n) = H{x,y{n)); 
Then for all y e such that k = J^Zi Vi, [H{x{t),y)] 
(C) Schrddinger equation with Bosonic Hamiltonian: If v : 



u{t]n{y)). 



X (Z 



solves: 



[!) For all n G {'L>of and t G M+, 



d , ^, 
-v{t-n) 



Hv{t;n), H={l-q) 



Vj + {l-q ^) ^ Sni=njq^ 'ffln^Vi 



(2) For all permutations of indices a G Sk, v{t; an) = v{t; ft); 

(3) For all n G (Z>o)'^ suc/i i/iai = 0, n) = for all t G 



(4) For all n G W^q^ ^{0; 



n 



Then for all y G Y'^ such that k 



H{x,y{n)); 



v{t;n{y)). 



Remark 2.8. The existence and uniqueness of global solutions to (A) is explained in Remark 12. 5i 
This then implies the existence of solutions in (C). It is not clear, a priori, that there exist solutions 
to (B). As we see in the proof of (B), the combination of the four conditions in (B) implies that 
restricted to n G W^q, u{t;n) = h{t;y{n)) for all t G M-(-. However, it is not clear that there exists 
a suitable extension of u outside the physical region W^q which satisfies the four conditions. Note 
that (B) should be considered as an advanced version of the method of images. Finally, though 
the above results are written for deterministic x (i.e., deterministic initial data) by linearity one 
can average over random x and achieve the same stated results with E^' [H{x{t), y)] replaced by its 
average over x, written as E,[H{x{t),y)], and the initial data for the ODEs likewise replaced by 
E[H{x,y)]. 

Proof. Call the three conditions contained in Definition 12.41 (A. 1) . (A. 2) and (A. 3). Part (A) follows 
from Theorem 12.21 since it implies that 



jE^[H{x{t),^. 



-TAZRPtox 



which matches (A.l). Along with this, the value of E^ [H{x{t),y)] is uniquely characterized by the 
initial data and the fact that (due to the definition of H) W [H{x{t),y)] = for all y G Y^ with 
yo > 0. Conditions (A. 3) and (A. 2) respectively match these properties and hence (A) follows. 

Part (B) follows by showing that if the four conditions for u given in (B) hold, then it implies that 
u{t; fi{y)) satisfies part (A) and hence that u{t] n[y)) = h{t; y). Thus we must show that (B) implies 
(A). Going between y and n notation, the initial data (A. 3) and (B.4) match, as do the conditions 
(A. 2) and (B.3). To check the system of ODEs (A.l), recall that the size of the cluster of elements 
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of n equal to i equals yj. Consider the cluster of elements equal to A^: ni = n2 = • • • = (every 
other cluster works similarly). In order to prove (A.l) it suffices to show that 

Vn 

(1 - g) ^ aArViu(t; ft) = apfil - qy^)Vy^u{t; n). (6) 



=1 

,A^,7V-l^ 



This is because Vyj^u{t;n) = u{t;fi{y ' )) — u{t; fi{y)) . Summing these terms over all clusters 
yields L''~'^^^^'^u{t-n{y)) and hence (A.I) follows. But (B.2) implies Viu{t;n) = qV^-'-V y^u{t;n) 
for i = I, . . . , T/TVi which implies ([6|). 

Part (C) also follows by showing that the combination of the four conditions for v imply that 
v{t\n[y)) satisfies (A) and hence v{t;n{y)) = h{t;y). As in (B), the initial data (A. 3) and (C.4) 
match, as do the conditions (A. 2) and (C.3). Also as in (B), it suffices to consider the cluster of A^. 
The portion of the Hamiltonian H corresponding to this cluster is 



Vn Vn 



Vn 



{l-q)aNY,q^''-'Vi, 



1=1 i<j 

where the equality follows from summing the factors involving each Vj. Due to the symmetry 
(C.2), Vjf (t; n) = S/yj^v{t; n) for alH G {I, . . . , yTv}. Hence the sum in i can be performed, yielding 
a7v(I — 1^'^)^yN- This is the same as in ([6|) and hence (C) follows as well. □ 

2.3. Nested contour ansatz solution. It is not a priori clear how one might explicitly solve the 
systems of ODEs in Proposition 12.71 Presently we show how this can be done for two distinguished 
types of initial data. 

Definition 2.9. For (7-TASEP, step initial data corresponds with ^^(O) = —i for i>l. 

For a G [0, 1) we say a random variable X is g-Geometric distributed with parameter a (written 
X ~ qGeo{a)) if 



W{X = k) = (a;, 



where (o; q)n = (I — a)(I — aq){l — aq^) • • • (I — aq"-~^) and (a; q)^ = (I — a)(l — aq){l — aq^) ■ ■ ■ . 
Half stationary initial data for g-TASEP corresponds with random initial locations for particles Xi 
for i > I given as follows: Let Xi ~ qGeo{a/ai) for i > I be independent; then set xi{Q) = —1 — Xi 
and, for i > 1, Xi = Xi^i — 1 — Xi. The result is that the gaps between consecutive particles i and 
i + 1 are distributed as (7-Geometric with parameter a/oj and are independent. When q = the 
step initial data is recovered (regardless of the Oj). 

Remark 2.10. When Oj = 1, the translation invariant measure on particle configurations in Z with 
independent gGeo(a) distributed distances between neighbors is an invariant or stationary measur^ 
for g-TASEP, cf. [8]. This explains the usage of the term half stationary (and likewise for ASEP). 

Theorem 2.11. Fix q G (0, 1), aj > for i > 1 and let n = (ni, . . . ,nfc). The system of ODEs 
given in Proposition \2. 7| (B) is solved by the following formulas. 
(I) For step initial data, 



. fc(fc-i) , , k 



l<A<B<k ^ j=l \m=l 

where the integration contour for za contains {qzB}B>A o-nd all Om's but not 0. 



Within the probabihty hterature, the term equilibrium is sometimes also used to describe such a measure, though to 
avoid confusion with the physical means of equilibrium statistical mechanics, we avoid this term. 
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(2) For half stationary initial data with parameter a > (such that aq ^ < am for all 1 < m < 
maxi{ni)), 

, fc(fc-l) k / i^j \ 

(2^0^ J J - qzs ^ i am - z, ) z, -a/q 

where the integration contour for za contains {qzB}B>A and all Om's but not a/q. 

On account of Proposition 12.71 and the uniqueness of solutions restricted tons W>q (see Remark 
12. Sp . the above formulas when restricted to n G ^>o imniediately yield the following (we will only 
state it for step initial data, though a similar statement holds for half stationary): 

Corollary 2.12. For q-TASEP with step initial data and n G W^q, 

k 



E 



. k{k—l) k / "^j \ 



where the integration contour for za contains {qzB^ b>a and all am's but not 0. 

Proof of Theorem \ 2. 1 1[ We need to prove that u{t; n) as defined in ([7]) satisfies the four conditions 
of Proposition [231(B). 

Condition (B.l) is satisfied by linearity and the fact that 

"|-(i-,K,v.]((n^).<-"=).o 

Condition (B.2) relies on the Vandermonde-like factors as well as the nested choice of contours. 
Without loss of generality assume that ui = n2- We wish to show that 

[Vi-qV2]u{t;n) =0. (10) 

Applying Vi — qV2 to the integrand in ([7]) brings down a factor of —a~^(zi — qz2). We must show 
that the integral of this new integrand is zero. This new factor cancels the denominator (zi — qz2) 
in 

ZA - ZB 



n (") 



On account of this we may deform the contours for zi and Z2 to be the same without encountering 
any poles. The term zi — Z2 in the numerator of (jlip remains and hence we can write 

u{t;fi)= / {zi - Z2)G{zi)G{z2)dzidz: 



2 

where G{z) involves the integrals m zs, . . . , Zk- Since the two contours are identical, this integral is 
clearly zero, proving (B.2). 

Condition (B.3) follows from simple residue calculus. When = there are no poles in the Zk 
integral at {am}m=i- Therefore, by Cauchy's theorem the integral is zero. 

Condition (B.4) likewise follows from residue calculus. Let us first consider the step initial data 
case of Theorem 12.111 This corresponds to initial data in (B.4) given by 

n(0; n) = H{x; y{n)) = 1 

for all n E W>o- Now consider ([8]) with t = 0. The zi contour can be expanded to infinity. The only 
pole encountered is at zi = {zi = cxd is not a pole because of the decay coming from am/ {am — Zj)). 
Because we pass it from the outside, the contribution of the residue is — times the same 
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integral but with every term involving zi removed. Repeating this procedure for Z2 leads to — 

and so on. Therefore the integral can be evaluated and canceling terms we are left with it equal to 

1 exactly as desired. 

Now consider the half stationary initial data case of Theorem l2.11[ We claim that this corresponds 
to initial data in (B.4) given by 

u{0;n)=E[H{x;y{n))]=Yl f[ / i " (l^) 

i=lm=n,+i + li=l 

Showing this requires a calculation. 

Lemma 2.13. Fix r > 1. If X is a q-Geometric random variable with parameter a G [0, 1) then 

1=1 ' 

SO long as aq~^ < 1; and otherwise the expectation is infinite. 

Proof. Using the g-Binomial theorem (see Section 17. ip , we may calculate 

which after canceling terms is exactly as desired. □ 

Recall that under half stationary initial data, the locations {xj(0)} are defined in terms of q- 
Geometric random variables {Xj}. Using this we have 

-Q qx„, (0)+n, =Y[q- S-=i =Y[ n ^'^^"^ • 

i=l i=l 1=1 m=ni+i+l 

Since the X's are independent, we can evaluate individually the expectation of g"*^™ using the 
above lemma, and we immediately find the right-hand side of (|12p . 

Now consider ([7|) with t = 0. As in the step initial data case, we successively peel off the contours 
and evaluate the effect via residue calculus. When we expand zi to infinity we now only encounter 
a pole at zi = a/q (which becomes zero when a = recovering the step initial data). Evaluating 
this residue we find 

ni , .fc„i (k-l)(k-2) k / rij \ J 

iii «™ - "/^ (2vr0'^-i J J J}^^^ ZA - qzs {^^^ a. - ) - a/q^ 

Expanding Z2 the pole is now at Z2 = Oi/q^ and a similar formula results from evaluating the residue. 
Repeating this procedure shows that n(0; n) is given by the right-hand side of (fT2]) as desired. □ 



3. A GENERAL SCHEME FROM NESTED CONTOUR INTEGRALS TO FrEDHOLM DETERMINANTS 

The output of Theorem 12.111 and Corollary 12.121 is that for step and half stationary initial data 
we have relatively simple formulas for a large class of expectations. In principle these expecta- 
tions should characterize the joint distribution of the locations of any fixed collection of particles 
Xni{t), . . . ,Xni{t) in g'-TASEP. One may hope to achieve this via certain generating functions. How- 
ever, the challenge is to find expressions for these generating functions which have clear asymptotic 
limits (in time and particle labels). For this, we focus here only on the distribution of a single 
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particle Xn (t) . Applying Corollary 12.121 with m = n yields a nested contour integral formula for 

There are two ways to deform this type of nested contour integrals so that all contours coincide. 
After accounting for the residues encountered during these deformations, we are led to two types 
of formulas for expectations: those involving partition-indexed sums of contour integrals and those 
involving single row-indexed sums of contour integrals. By taking suitable generating functions of 
these indexed sums of contour integrals we are led to two types of Fredholm determinants. All 
of these manipulation are quite general and can be done purely formally. Given some analytic 
estimates these manipulations turn into numerical equalities as is the case for g-TASEP. 

We record these manipulations (and conditions for them to hold as numerical equalities) as well 
as their consequences without proofs, since they can be found in Section 3.2 of [8]. We do this for 
completeness and also because when we turn to consider ASEP, the same manipulations will be 
used. When we apply this to g-TASEP we will only consider step initial data. If we consider half 
stationary for any a > fixed, then when k gets so large that a > q^, the expectation E 
is infinite. Thus, when forming a generating function from these g-moments, we are forced to take 
a = 0, which corresponds to the step initial data. 

Before going into these manipulations the reader may want to quickly browse Section 17.11 where 
we record some useful g-deformations as well as briefly review Predholm determinants. 

3.1. Mellin-Barnes type determinants. The below proposition describes the result of deform- 
ing the contours of a general nested contour integral formula in such a way that all of the poles 
corresponding to za = qzs for A < B are encountered. The residues associated with these poles 
group into clusters and hence the resulting formula is naturally indexed by partitions A = (Ai > 
^2 > • • • > 0). Notationally, we write A h /c if Aj = A;, A = l'"i2'"2 ... if ^ appears rrii times in A 
(for all i > 1), and ^(A) = for the number of non-zero elements of A. 

Definition 3.1. For a meromorphic function f{z) and k > 1 set A to be a fixed set of poles of / 
(not including 0) and assume that q^K is disjoint from A for all m > 1. Define 

where the integration contour for za contains {qzB}B>A, the fixed set of poles A of f{z) but not 
or any other poles. 



Proposition 3.2. We have that for /z^ as in Definition \3.1 



ei\) ^(A) 



1 (I — a) f f 1 ^ ' 

f^k = k,\ Yl , , T^r^ / • • • / det n fi^^)f(l^j) ■ ■ ■ fiQ'^-'wj)dw, 

A=1™12™2... 

(14) 

where the integration contour for wj contains the same fixed set of poles A of f and no other poles. 
Proof. This is proved in [8] as Proposition 3.2.1 via residue calculus. □ 

As a quick example consider f{z) which has a pole at z = 1. Then the z^-contour is a small 
circle around 1, the z^-i-contour goes around 1 and q, and so on until the zi-contour encircles 
{l,q, . . . ,q^~^} (this is illustrated for A; = 3 in Figure [2]). All the w contours are small circles 
around 1 and can be chosen to be the same. 

We form a generating function of the /x^ and identify the result as a Fredholm determinant. 
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Figure 2. Possible contours when k = 3 for the zj contour integrals in Proposition 13.21 



Proposition 3.3. Consider fi^ as in equation defined with respect to the same set of poles A 
of f{w) for k = 1,2,... and set Ca to be a closed contour which contains A and no other poles of 
f{w)/w. Then the following formal equality holds 

k>0 ^' 

where det(/ + J^^J) is the formal Fredholm determinant expansion of : L^(Z>o x Ca) — )• L^(Z>o x 
Ca) defined in terms of its integral kernel 

(1 - qrK'''fim)fiqwi) ■ ■ ■ f{q^'-^wi) 



q^^Wi — W2 



The above identity is formal, but also holds numerically if the following is true: For all w, w' £ Ca 
and n > 1, \q'^w — w'\~^ is uniformly bounded from zero; and there exists a positive constant M 
such that for all w £ Ca and all n>0, |/(g"w)| < M and |(1 - q)C\ < M'^ 

Proof. This is proved in [8], Proposition 3.2.8. The proof amounts to reordering the sums defining 
fik and recognizing a Fredholm determinant. □ 

We may replace the space -L^(Z>o x Ca) by L?'{Ca) via the following Mellin-Barnes representation. 

Lemma 3.4. For all functions f which satisfy the conditions below, we have the identity that for 

Ce{C:|CI<i,C^K+}. 



oo 

E = / n-sw + s){-crfiqids, 



(15) 



where the infinite contour Ci^2,... is a negatively oriented contour which encloses 1,2,... and no poles 
of /(<?*) (e.g. Ci^2,... = ^ + iM oriented from ^ — loo to ^ + loo), and z^ is defined with respect to 
a branch cut along z £ M". For the above equality to be valid the left-hand-side must converge, and 
the right-hand-side integral must be able to be approximated by integrals over a sequence of finite 
contours which enclose the poles at 1,2, ... ,k and which partly coincide with Ci^2,... such a 
way that the integral along the symmetric difference of the contours Ci^2,... and Ck goes to zero as 
k goes to infinity. 
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Figure 3. Left: The contour Dji^^.j^; Right: The contour D^^- 
Proof. The identity follows from Resr(— s)r(l + s) = (—1)'^"'"-'^. □ 

s=k 

Definition 3.5. The infinite contour Djid is defined as follows. Dji^^ goes by straight lines from 
R — LOO, to R — id, to 1/2 — id, to 1/2 + id, to R + id, to R + ico. See Figure [3] for an illustration. The 
finite contour Dji ^-k is defined as follows. Let p,p be the points (let Im(p) > 0) at which the circle 
of radius k + 1/2, centered at 0, intersects Dr ^. Then D^ ^-k is the union of the portion of -D^.d 
inside the circle with reversed orientation, with the arc from p to p (oriented counter-clockwise). 

Proposition 3.6. Assume f{w) = g{w) / g{qw) for some function g. Then the following formal 
equality holds 

det{I + Kl) = det{I + K^) 

where det(I + K^) is given in Proposition \3.3\ and where det(/ + is the formal Fredholm de- 
terminant expansion of : L'^{Ca) L'^{CjC}. The operator K'^ is defined in terms of its integral 
kernel 

Kl{w,w') = ^ [ r(-.)r(i + .)(-(i-g)c)^ , ^ , ds. 

2Tri Jci,2,... 9{q w) q w - w' 
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The above identity holds numerically if det{I + K^) is a convergent Fredholm expansion and if 
Ci,2,... is chosen as Dfi ^ with d > and R > such that 



inf \q^w — w'\ > 0, and sup 

w,w'eCA w,w'eCA 

In that case the function ^ i— )■ det(/ + K?) is analytic for all ^ ^ IR+. 



< oo. 



Proof. This result can readily be extracted from the proof of [8] Theorem 3.2.11. In fact, the strong 
analytic bounds which we require can be significantly relaxed, however, as they will be sufficient for 
our purposes, we do not explore this. □ 

We say that Fredholm determinants similar to det(J + K'^) are of Mellin-Barnes type (because 
(jlSp is a basic tool for classical Mellin-Barnes integrals). 



3.2. Cauchy type determinants. Instead of deforming contours so as to encounter the za = QZb 
poles, we may deform our contours to only encounter the pole at 0. The residue calculus becomes 
easier and the resulting sum of contour integrals is indexed by partitions with just a single row 
(equivalently by nonnegative integers). 

Definition 3.7. For a meromorphic function f{z) and k > 1 set A to be a fixed set of poles of / 
and assume that q"^A is disjoint from A for all m > 1. Define 

{2irt)^ J J ZA — QZb Zi 

^ ' •> •> l<A<B<k ^ ^ ^ i=l * 



where the integration contour for za contains {(IZb} b>a^ the fixed set of poles A of f{z) and 0, but 
no other poles. 

Notice that and /i^ differ only by the inclusion of in the contour for fi^. They can be related 
via the following. 

Proposition 3.8. Assume /(O) = 1. Then, 

Proof. This is proved in [H], Proposition 3.2.5. □ 

If, for instance, we assume now that A contains all poles of /, then we can deform the contours 
in (I16p to all lie on a single, large circle. The following symmetrization proposition then applies. 



Proposition 3.9. If the contours of integration in / fi6]) can he deformed (without passing any poles) 
to all coincide with a contour Ca, then 



IJ-k 



fc,! (l-r^)^ r r . .k k 

k\ (27r6)'= 



/ •••/ det 3^ f\f{wj)dwj. (17) 

Proof. This is proved in [8], Proposition 3.2.2. □ 
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Proposition 3.10. // the contours of integration in U6\} can be deformed (without passing any 
poles) to all coincide with a contour Ca, then the following formal identity holds 

^/ifcf- = det(/ + Ci<'') 

Kq. 



k>0 



where det{I + K^) is the formal Fredholm determinant expansion of : L'^{Ca) L'^{Ca) defined 
in terms of its integral kernel 



qw' — w 



The above identity is formal, but also holds numerically for ( such that the left-hand side converges 
absolutely and the right-hand side operator is trace-class. 

Proof. This is proved in [8], Proposition 3.2.9. □ 

Remark 3.11. By considering the Fredholm series expansion (whose terms are given by (|17p ) it 
is clear that since / arises multiplicatively, it can be paired either with Wi or Wj in the Cauchy 
determinant. As a consequence, it follows that 

det(/ + QK^) = det(/ + QK'^) 

where : LP'{Ca) — t- L'^{Ca) is defined in terms of its integral kernel 

. fH 



K\w,w') = (1 



qw — w' 



We call Fredholm determinants of this form Cauchy type. 

3.3. Application to g-TASEP. The following theorems about g-TASEP are applications of the 
manipulations of the previous section. The required estimates necessary to make these numerical 
equalities are provided in [8]. 

3.3.1. Mellin-Barnes type Fredholm determinant for q-TASEP. 

Theorem 3.12. Fix Q < q <1 and n > 1. Fix < (5 < 1 and ai, . . . , a„ such that for all i, Oi > 

and \ai — 1| < d for some constant d < Consider q-TASEP with step initial data and jump 

parameters ai. Then for all t G R+ and G C \ the following characterizes the distribution of 
Xn{t): 

^ + Kr^^''''') (18) 



E 



{Qq-^it)-q)^ 



where dei{I -\- K'^ tasep-^ Fredholm determinant of : L'^{Ca) — L'^{Ca) for Ca a positively 

oriented circle \w — 1\ = d. The operator is defined in terms of its integral kernel 



1 pLOO + O 

K,-TASEP^^^^,^^ r(-s)r(i + s)(-g-"c) 

^ 27ri gil'w) q'w - w' 



ds 



where 



9{w) {w/am;q)oo'^ ' 

m=l 

Proof. This is proved in [8], Theorem 3.2.11. A similar approach is described in its entirety in the 
proof of Theorem 15.31 for ASEP. 

□ 
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The above is an e^-Laplace transform and can be inverted via Proposition 17.11 Since is 
supported on {— n, — n + 1, . . .}, in order to apply Proposition 17. 1 1 it is necessary to shift everything 
by n. Let /'^(C) = det(/ + /^9~'^^sep^ redefine Cm to encircle the poles C = Q~^^ for —n < 
M < m — n. Under these modifications, Proposition 17.11 gives P(xn(t) = m). 

3.3.2. Cauchy type Fredholm determinant for q-TASEP. 

Theorem 3.13. Fix 0<q<l,n>l and ai, . . . ,an such that for all i, Oj > 0. Consider q-TASEP 
with step initial data and jump parameters > for all i > 1. Let Xn{t) by the location of particle 
n at time t. Then for all C & C 

det(J + Ci^^-^^^^0 

(C;^)oo ^'"^ 



E 



1 



[(q^nit)+n.q^^ 



where det{I + (K'^ TASEP^ entire function of C, and is the Fredholm determinant of K'^ TASEP . 

L'^{Ca) — ^ L'^{Ca) defined in terms of its integral kernel 



K1-TASEP^^^^>^_ f{w 



qw' — w 



with 



and Ca a star-shaped contour with respect to (i.e. it strictly contains and every ray from 
crosses Ca exactly once) contour containing ai, . . . , a„. 

Proof. This is proved in [8], Theorem 3.2.16. A similar approach is described in its entirety in the 
proof of Theorem [531 for ASEP. □ 

The above shows that det(/ + C,K^~'^^'^^^) / q)oo equals the e^-Laplace transform of 

4. Duality and the nested contour integral ansatz for ASEP 

The asymmetric simple exclusion process (ASEP) was introduced by Spitzer [36j in 1970 and also 
arose in biology in the work of MacDonald, Gibbs and Pipkin [26j in 1968. Since then it has become 
a central object of study in interacting particle systems and non-equilibrium statistical mechanics. 

The ASEP is a continuous time Markov process with state r]{t) = {rix{t)}xez £ {0, 1}^ at time 
t >0. The r]x{t) are called occupation variables and can be thought of as the indicator function for 
the event that a particle is at site x at time t. The dynamics of this process is specified by non- 
negative real numbers p < q (normalized by p + g = 1) and uniformly bounded (from infinity and 
zero) rate parameters {a^jxez- For each pair of neighboring sites {y, y + 1) the following exchanges 
happen in continuous time: 

7] r/^'S^+i at rate Uyp if {i]y,r]y+i) = (1,0), 

■q h-^ T^^^'f+i at rate ayq if {r]y,r]y+i) = (0, 1), 

where jf^y^^ denotes the state in which the value of the occupation variables at site y and y + 1 
are switched, and all other variables remain unchanged. All exchanges occur independently of each 
other according to exponential clocks. These dynamics are called the ASEP occupation process and 
are defined in terms of the generator U"^'^ which acts on local functions / : {0, 1}^ — t- M by 

{L°'"'f){ri) = ^ay [pr^yil - r^y+,) + q{l - Vy)vy+i] (/(^?^'^+^) - /(r/)) . (21) 
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Figure 4. ASEP with four particles: xi = 5, X2 = 6, 2:3 = 9, and X4 = 11. The first 
two particles form a cluster, and the third and fourth form two separate clusters. 
The arrows represent admissible moves. 

The existence of a Markov process this generator is shown, for example, in |24| . 

In terms of particles, the dynamics of ASEP are that each particle attempts, in continuous time, 
to jump right at rate pay and to the left at rate goy-i (presently the particle is at position y G Z), 
subject to the exclusion rule that says that jumps are suppressed if the destination site is occupied. 
We assume p < q (drift to the left) and define 7 := g — p > and t := p/q < 1. 

The ASEP preserves the number of particles, thus we can consider ASEP with k particles as a 
process on the particle locations. Define 

W'' = {x = {xi,X2, . . . , Xfc) e : Xi < X2 < . . . < Xk} 

and xf = {xi,. . . ,Xi-i,Xi ± l,Xi+i, . . . ,Xk). Then x(t) = (xi(t) < X2{t) < ■ ■ ■ < Xk{t)) G 
denotes the locations of the k particles of ASEP at time t. 

In order to describe the generator of ASEP in terms of particle locations it is convenient to intro- 
duce particle cluster notation. A cluster is a collection of particles next to each other: xi — — 
1 = • • • = Xj+j — j. There is a unique way of dividing the particles x into clusters so that 
each cluster is separated by a buffer of at least one site: Let c(x) be the number of such clus- 
ters, i{x) = {ii, . . . be the collection of labels of the left-most particles of each cluster, and 
r(x) = (ri, . . . ,rc) be the collection of labels for the right-most particles of each cluster. For in- 
stance, if /c = 4 and xi = 5, X2 = 6, X3 = 9, x^ = 11 then c{x) = 3, i{x) = (1,3,4) and 
r(x) = (2,3,4). _ 

For A; > 1, the ASEP particle process generator acts on bounded functions / : — )• M by 

(LP-V)(^)= E a,.-iP[f{K)-f{x)]+ a,,g[/(f+)-/(f)]. 

We will consider initial configurations for ASEP in which there is at most a finite number of non- 
zero occupation variables (i.e. particles) to the left of the origin - we call these left-finite initial data. 
When ASEP is initialized with left-finite initial data, its state remains left-finite for all time (simply 
because it will always have a left-most particle). We will use E'' and P'' to denote expectation and 
probability (respectively) of the Markov dynamics on occupation variables with initial data r] (and 
likewise and for the Markov evolution on particle locations with initial data x). When the 
initial data is itself random, we write E and P to denote expectation and probability (respectively) 
of the Markov dynamics as well as the initial data. We also use E and P when the initial data is 
otherwise specified. 

4.1. Duality. The following result shows that ASEP with general bond rate parameters is self-dual. 
Recall that t = p/q < 1 by assumption. Define the following functions of a state 77: 

Nxit) = V ^lyit), Qx{t) = T^^W, Qxit) = Q^'W-Q^-i(^) = r^^-Wr?.(t). (22) 

^ — ' T — 1 

y=-oo 

Theorem 4.1. Fix non-negative real numbers p < q (normalized hy p + q = \) and uniformly 
hounded (from infinity and zero) bond rate parameters {ax}xez- For any k > 1, the ASEP occupation 
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process 7]{t) with state space {0,1}^, and the ASEP particle process x{t) with state space are 
dual with respect to 

k 

H{v,S) = Y{Q.^{t). 

i=l 

If we restrict to ax = 1 we can prove another ASEP self-duality. 

Theorem 4.2. Fix non-negative real numbers p < q (normalized by p + q = 1) and bond rate 
parameters Ux = 1- For any k > 1, the ASEP occupation process ri{t) with state space {0, 1}^, and 
the ASEP particle process x{t) with state space are dual with respect to 

k 

H(n,x) = l[QxM- 
1=1 

Recall that the concept of duality is given in Definition 12.11 A few remarks are in order. 

Remark 4.3. For p < q, both forms of duality are trivial for initial data which is not left-finite, 
since then Qxit) = and likewise Qx{t) = 0. 

Remark 4.4. For the symmetric simple exclusion process {p = q = 1/2) the duality from Theorem 
14.11 has been known for some time (see [21] Chapter 8, Theorem 1.1). For p < q, and under the 
assumption that all Ux = 1 (i.e. no bond dependence in jump rate parameters) the result of Theorem 
14.11 was discovered by Schiitz [35] in the late 1990's via a spin chain representation of ASEP. The 
approach used therein to show duality was computationally based on a Uq{sl2) symmetry for ASEP 
[31] . Our results apply for general bond rate parameters and the proof proceeds directly via the 
Markov dynamics, without any use of, or reference to, the Uq{sl2)- 

The duality of Theorem 14.21 appears to be new. It does not seem possible to extend it to general 
ax- For instance, when A; = 1, as a function of the process 7?(t), H{ri{t),x) only changes value when 
a particle moves across the bond between x and x + 1. This only involves the rate Ox- On the other 
hand, as a process of x{t), H{rj,x(t)) changes value when a particle moves across either the bond 
between x — 1 and x, or the bond between x and x -|- 1. This involves the rates a^-i and a^. Hence 
the two sides can only match when ax-i = Ox- 

Remark 4.5. Gartner [16] observed that ASEP respected a microscopic (i.e., particle- level) version 
of the Hopf-Cole transform (see e.g. the review [l3]). This observation is equivalent to the /c = 1, 
ax = 1 case of the duality given in Theorem 14.21 It says that 

dQx{t) = {pQx-i{t) + qQx+i{t) - Qx{t)) dt + QxdM{t) 

where dM{t) is an explicit martingale. This is a particular semi-discrete SHE (different than the 
one coming from g-TASEP, Definition 16. ip with a somewhat involved noise (the martingale is not 
exactly a discrete space-time white noise). A Feynman-Kac representation for this equation shows 
that Qxit) can be thought of as a polymer partition function with respect to an environment defined 
by the martingale. Therefore, Theorem 14.21 can be thought of as a version of the polymer replica 
approach (see Section [6?2]) . 

The proof of the two duality theorems boils down to two propositions which we now state and 
prove. After this we prove the theorems. 

Proposition 4.6. Fix non-negative real numbers p < q (normalized by p + q = \) and uniformly 
bounded (from infinity and zero) bond rate parameters {ax}x&- Then, with rj, x, and H{rj, x) defined 
in Theorem \4.1\ 

L'"''H{r], x) = LP'^'-'Hirj, x), (23) 
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where the generator L"^'^ acts in the rj variable and the generator acts in the x variable. 



Proof. We will first prove the desired result for a single cluster configuration x = + + 
and then easily deduce it for general x E W^. For the single cluster x, by the definition of L"'^^, 

i=-l 

where 

We now compute the Ay^s explicitly. There are three different types of Ay that must be considered: 
(1) Aa;-i; (2) Ax+i for < i < £ - 1; (3) A^+i,. We use the notations (l22]l . 
(1) Consider A^-i. We may rewrite 

H{i], x) = r^--V^--i??, J] Q,+„ and ^(r?-i'-, f) = r^^" V"-??,-! J] Q,+j. 

j=i i=i 

Thus 



i=i 

(2) Consider Ax+i for < i < £ — 1. We may rewrite 



(24) 



H{i], x) = r 



Vx+iVx+i+l Qx+j) 
j=0 



and 



Thus 



i=o 



■x+t 



n ^- 



\jyi,i+i / 



X {prix+i{l - rjx+i+i) + g(l - ??x+i)^x'+j+i) [r''"+'+i - r^'^+'jry^+jry^+j+i. 
(3) Consider A^+i. We may rewrite 



-1 



-1 



,X 



T^'^+'-'Vx+l+lllQ. 



(25) 
(26) 



x+j- 



j=0 



j=0 



Thus 



£-1 



\3=0 

X (pr/a;+£(l - 77x+£+i) + q{l - r]x+£)r]x+£+i) [nx+e+i - Vx+e] ■ 



(27) 
(28) 
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Observe that A^+i = for < i < £ — 1. To see this, it suffices to consider the four values that 
the pair (?7a;+i, ?7^.+i+i) may take: For (0,0) or (1, 1) 

{pr]x+i{l - Vx+i+i) + q{l - 7]x+i)Vx+i+i) = r''^+' - r''-+'+i = 

and thus (125]) = 0; for (0, 1) or (1, 0), the factor r]x+ir]x+i+i = and thus (f25]) = again. 
The above observation shows that, in fact, 

[1°'"' f){it]) = ax-iAx-i + Qx+iAx+t. 

In hght of equations (|24p and (j27p we may rewrite 

Ax-i{i]) = M(7?X._i(77) and Ax+d^) = M{r^)A',^,{r,), 

where 
and 

^'x-iijl) = r]xu{pVx-i{l-r]x)+q{l-Vx-i)Vx)[r'^''Vx-i-T'^''-^r]x] 

A'^+i = t''--'??x {prix+eC^ - r]x+e+i) + ^(1 - Vx+e)Vx+e+i) [rix+£+i - Vx+e] ■ 

Now turn to the right-hand side of equation (|23p . We may also factor M{rj) out from that 
expression 

RHS(I23]) = M(ry) [ax-iprjx-iVx+e + ax+eqT'^''-^r]xT'^^+''r]x+e+i - {ax-iQ + ax+^p)T"''"~'^x??x+f] • 
Therefore, for the single cluster case of the proposition, we are left to prove 

ax-iA'^_.^ + ax+eA'^+^ = ax-iprjx-irjx+e + ax+£gT"''""'??xT"''"+'??x+^+i - {ax-iQ + ax+^p)-^''""'??x^x+£- 

The above equation is a function of only four occupation variables ry^-ii ^a;+£ and r]x+e+i and 
one can systematically check that for all sixteen combinations of values of these variables, the above 
equation is true. In fact, it is even easier than this since the coefficients of a^-i and ax+i coincide 
separately. For instance we must show that A'^_i = rjx+i {prjx-i — qr^'^^^rjx). There are only four 
cases of {rjx-i,rjx) we have to be considered and this can be confirmed in one's head (similarly for 
^x+f)- This completes the proof of Proposition 14.61 for x with just a single cluster. 

From a general x S W'' there may be many clusters, each pair separated by at least one empty 
site. The terms in H{rj, x) factor into clusters and the generator L°^'^ acts on each of these clusters 
according to the above proved single cluster result. This immediately yields the general statement 
and completes the proof. □ 

Proposition 4.7. Fix non-negative real numbers p < q (normalized by p + q = 1) and set all bond 
rate parameters = 1. Then, with r], x, and H{r],x) defined in Theorem \4-S\ 

L"''H{r,,x) = LP'''~'H{r],x), 

where the generator L"'^^ acts in the rj variable and the generator L^""^^ acts in the x variable. 

Proof. As in the proof of Proposition 14.61 we will first prove the desired result for a single cluster 
configuration x = {x,x + I, . . . ,x + £) and then easily deduce it for general x G W^. For the single 
cluster X, by the definition of L°^'^, 

I 

L°^'H{r],x) = Y,Ax+i{v), 

1=0 
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where 

MV) = {Privi^ - Vy+i)+q{l - r]y)Vy+i) [H{r,y^y+\x) - H{r],: 
By groupmg terms this may be rewritten as 



Ax+iiv) = Qx+i~l 



\ 



, j=0 , 



j=0 



j=0 



In order to get the second hne above, we utihzed the definition of Qx+i-i and separately the fact 
(which can readily be checked) that for the four possible pairs of values that (^Tjx-\- i,r]x+i+i) can take, 

{pVx+i{l - Vx+i+i) + qVx+i+i{l - r]x+i)) [r^-+^+i - t''-+'] =p + qT^-+- t'^-+^+^ - t'^-+\ 

Recall that we seek to show 

Ax+i = pQx~i JJ Qx+j + qQx+£+l llQx+j-llQx+j. 

i=0 j=l j=0 j=0 



Factoring out (Jlfco Qx+j-i ) from both sides we are left to prove 



^+1 



(p _j_ qq-Vx + i ^^x + i+l jVx + i^ — p J~J q-Vx+j _j_ q J~J q-Vx+j J~J q-Vx+j _ 

j=l j=o j=o 



(29) 



The terms in the left-hand side of the above expression can be grouped as 



i=l j=0 



j=0 



j=0 



-r]x+j 



We may now utilize the easily checked identity that for r/ G {0, 1}, 

p + gT2^-r'' = 0, 

to see that the above expression reduces to the right-hand side of ([29|) . thus completing the proof 
of Proposition 14.71 for x with just a single cluster. 

From a general x G W'^ there might be many clusters, each pair separated by at least one empty 
site. The terms in H{ri, x) factor into clusters and the generator L°'^'^ acts on each of these clusters 
according to the above proved single cluster result. This immediately yields the general statement 
and completes the proof. □ 

Before giving the proof of Theorem 14. II we define the following system of ODEs. 

Definition 4.8. We say that h{t\ x) : M-|_ x W'' — )■ M. solves the true evolution equation with initial 
data ho{x) if 

(1) For aU x G W'' and t G M+, 

^h(t;x) = LP''''h(t;xy, 
at 
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(2) There exist constants c, C > and 6 > such that for ah x G W^, t G [0, 6] 

\h{t;x)\ < Ce^ll^lli; 

(3) As t — 0, h{t;x) converges pointwise to ho{x). 

Proposition 4.9. Assume that there exists constants c, C > such that for all x G , 

\hQ{x)\<Ce^\\^\\\ (30) 

Then there exists a unique solution to the system of ODEs given in Definition \4.8\ which is given by 

:=E-*'^[/io(x(0))] (31) 

where the expectation is with respect to the ASEP particle process x(-) started at time —t in config- 
uration x. 

This existence and uniqueness result is proved in Section [8l We use this result presently in the 
proof of Theorem 14.11 and also later in the proof of Theorem 14.131 It is in the second application 
of this result that we fully utilize the weakness of conditions 2 and 3 in Definition 14.81 

Our present 



Proof of Theorem \4.1\ We follow the same approach as in the proof of Theorem 
theorem follows from Proposition 14.61 along with Proposition 14.91 Observe that 



d 
'dt 



L°==i/(r?(t),f) 
LP''''H{'n{t),x 



The equality of the first line is from the definition of the generator of rj{t); the equality between 
the first and second lines is from the commutativity of the generator with the Markov semi-group; 
the equality between the second and third lines is from applying Proposition 14.61 to the expression 
inside the expectation; the final equality is from the fact that the generator LP^''* now acts on the 
X coordinate and the expectation acts on the rj coordinate. This shows that, as a function of t 



and X, E'' 



H{r]{t),x) solves the true evolution equation of Definition 14.81 (checking condition 2 is 



straightforward and condition 3 can be checked as in the proof of Proposition 14.9P . 



On the other hand, Proposition 14.91 implies that E'' H{r],x{t)) 
equation of Definition 14.81 and that it is the unique such solution. ' 
to show the claimed duality. 

Proof of Theorem This follows exactly as in the proof of Theorem 14.11 with Proposition 14.61 
replaced by Proposition 14.71 □ 



also solves the true evolution 

'his proves the desired equality 

□ 



4.2. Systems of ODEs. As a result of duality, we provide two different systems of ODEs to 



characterize E^ 



H{rj{t),x) 
ODEs associated to g-TAS 



. These two systems should be compared to the first two systems of 
IP duality, given in Proposition 12.71 



We first state the result for the H{r], x) duality. 
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Proposition 4.10. Let rj be a left-finite occupation configuration in {0, 1}^ and r]{t) be ASEP 
started from r]{0) = rj. 

(A) True evolution equation: Ifh{t; x) : x — ^ M solves the system of ODEs given in Definition 
with initial data hQ{x) = H{r],x), then for all x E VF^, H{r]{t),x) = h{t;x). 



(B) Free evolution equation with boundary condition: If u : IR4. x — )■ M solves: 

(1) For all X and t € R+, 

d ^ 

—u{t;x) = ^ [a^^-lpu{t■,x~) + ax^qu{t;xf) - {a-j^^^iq + paxju{t; x)] ; (32) 

i=l 

(2) For all x ^ Z,^ such that for some i G {1, . . . ,k — 1}, Xj+i = Xj + 1, 

pu{t; xr^-^) + qu{t; xf) = u{t; x); (33) 

(3) There exist constants c, C > and 6 > such that for all x E W^, t G [0, 5] 

\u{t;x)\ < Ce^ll^lli; 

(4) For all x G W'', as t ^ 0, u{t;x) H(j],x). 
Then for all X <^W^,Wi H{r]{t),x) =u{t;x). 

Proof. Part (A) is an immediate consequence of the duality result of Theorem 14.11 along with its 
proof. Call the three conditions contained in Definition 14.81 (A.l). (A. 2) and (A. 3). 

Part (B) follows by showing that if the four conditions for u given in (B) hold, then it implies 
that u{t;x) restricted to x G actually satisfies conditions (A.l), (A. 2) and (A. 3). Conditions 
(B.3) and (B.4) immediately imply conditions (A. 2) and (A. 3), respectively. It is easy to check that 
the boundary condition (B.2) along with the free evolution equation (B.l) combine to yield the 
generator LP^''* and hence yield (A.l). Applying part (A) we see that given the conditions of (B), 

we may conclude that for all x G W^, E** H{rj{t),x) = u{t;x). □ 

We have an almost identical result and proof associated with the H{r], x) duality. 

Proposition 4.11. Let r] be a left-finite occupation configuration in {0, 1}^ and r]{t) be ASEP 
started from r/(0) = rj. 

(A) True evolution equation: Ifh{t; x) : R_). x — )■ R solves the system of ODEs given in Definition 
\4.8\ with initial data ho{x) = H{r], x), then for all x G W'' , W [H{ri{t),x)] = h{t; x). 

(B) Free evolution equation with boundary condition: If u : x — t- M solves: 

k 



(1) For all X £l} and t G 



^u{t; x) = ^ [pu{t; Xi ) + qu{t; f+ ) - u{t; x)] ; 
1=1 

(2) For all x Z,^ such that for some i G {1, . . . ,k — 1}, Xj_|_i = Xj + 1, 

pu{t; J + qu{t; xf) = u{t; x); 

(3) There exist constants c, C > and 6 > such that for all x G W'^ , t G [0, 6] 

Ht;x)\ < Ce^ll^lli; 
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(4) For all x G W'' , as t ^ 0, u{t;x) H{i],x). 
Then for all x e , [H{r]{t),x)] = u{t;x). 

Proof. Similar to that of Proposition 14.101 □ 

4.3. Nested contour integral ansatz. From now on we assume that all bond rate parameters 
Oa; = 1, in which case equation (j32|) becomes 

d ^ 

— x)=^ [pu{t; x~) + qu{t;xf) - u{t; x)] . (34) 

i=l 

It is not a priori clear how one might explicitly solve the systems of ODEs in Propositions 14.101 
and 14. Ill For g-TASEP, when confronted with the analogous problem of solving the system of ODEs 
in Proposition 12.71 we appealed to a nested contour integral ansatz which was suggested from the 
algebraic framework of Macdonald processes (into which g-TASEP fits). 

ASEP, on the other hand, is not known to fit into the Macdonald process framework, nor any 
similar framework from which solutions to these systems of ODEs would be suggested. Nevertheless, 
we demonstrate now that we may apply a nested contour integral ansatz. We focus on solving the 
system of ODEs in Proposition 14.101 for two distinguished types of initial data. Notice that in the 
below theorem, the contours are not nested, however they are chosen in a particular manner to 
avoid poles coming from the denominator za — tzb- 

Definition 4.12. For p G [0, 1] consider an i.i.d. collection {yE}a;>i of Bernoulli random variables 
taking value 1 with probability p. Then the step Bernoulli initial data for ASEP is given by setting 
Vxi^) = for X < and r]x{0) = 1^ for x > 1. When p = 1 this is called step initial data and 
(deterministically) rj^^O) = lx>i- We also define 6 = p/{\ — p). 
Define the function 



fx;p{z) = exp 



' (35) 



1 + z/t J T + Z Z 



{l + z)ip + qz) 

When p = 1 (and hence 6 = oo, the definition of fx;iiz) corresponds to the expression above, with 
the final fraction removed. Also define 



Fx;piz) = exp 



^2 



z{p - q) _^ 



l + z 
1 + z/t 



(36) 



{z + l){p + qz) 

and likewise extend to p = 1. 

Finally, define an integration contour C--r--i as a circle around — r, chosen with small enough 
radius so that —1 is not included, nor is the image of C„t-;-i under multiplication by r. It is also 
important that t6 and are not contained in C_t-;-1) but these facts are necessarily true from the 
definition. 

Theorem 4.13. Fix non-negative real numbers < p < q (normalized hyp + q = Ij and set all bond 
rate parameters = 1. Consider step Bernoulli initial data with density p G (0, 1]. The system of 
ODEs given in Proposition \4-10\ (B) is solved by the following formula 

fc(fc-i) ^ 

uit;x) = ■■■ n — —T\fxrA^i)dzi, (37) 

y I J \<A<B<k ^ ^ i=\ 

where the integration contour is given by C~t-\. 

As an immediate corollary of Theorem 14.131 and Proposition 14.101 (B) we find formulas for joint 
moments of the Qxi^^ defined in (j22]) . 
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Corollary 4.14. Fix k > 1, non-negative real numbers < p < q (normalized by p + q = 1) and 
set all bond rate parameters o^,- = 1. For step Bernoulli initial data with density p G (0, 1] and any 

k 



fc(fc-i) 



E 



Qx^it) ■ ■ ■ Q,,it) 



where the integration contour is given by C- 



n 



l<yl<B<A: 



ZA - ZB 
ZA - TZb 



Y{fx,;p{^i)dZi 



Proof of Theorem \4-13 . We give the proof for step Bernoulli initial data with p G (0, 1), and hence 
6 G (0, oo). The modification for the p = 1 case is trivial. 

We need to prove that u{t]x), as defined in (|37p . satisfies the four conditions of Proposition 14. lUl 
(B). 

Condition (B.l) is satisfied by linearity and the fact that 



dt 



fx;p{z) = 0, 



where AP''^5(x) = pg{x — 1) + qg{x + 1) — g{x) acts on the x- variable in fx-p{z). 

Condition (B.2) relies on the Vandermonde factors as well as the choice of contours. Without 
loss of generality assume that X2 = xi + 1. We wish to show that 

pu{t; X2) + qu{t\ x^) — u{t; x) = 0. 

Thinking of the left-hand side as an operator applied to u{t;x), we compute the effect of this 
operator on the integrand of (j37p and find that it just brings out an extra factor in the integrand 
(when compared to n(t;x^)) which is 

l + zi\fl + Z2\ f 1 + Z2\ ^ _^ {p-q)/t 



p + q 



(zi -TZ2)- 



(38) 



1 + Zi/t J \1 + Z2/T J \1 + Z2/tJ ' {1+zi/t){1 + Z2/t) 

We must show that the integral with this new factor times the integrand in (j37p is zero. The factor 
{zi — TZ2) cancels the term corresponding to j4 = 1 and B = 2 m. the denominator of 

ZA - ZB 



n 



ZA - TZB 



l<A<B<k 

The term {zi — Z2) in the numerator remains, and the additional terms coming from (j38p 
symmetric in zi and Z2- Therefore we can write 



are 



u{t; x) 



{zi - Z2)G{zi)G{z2)dzidz2 



where G{z) involves the integrals in 2:3, . . . , z^. Since the contours are identical, this integral is zero, 
proving (B.2). 

Condition (B.3) follows via very soft bounds. Observe that as z varies along the contour C_t;-i, 
and as t varies in [0, 5] for any 5, it is easy to bound \ fx:p{z)\ < Ce'^H^H^ for some constants c, C > 0. 
Since the contours are finite and since the other terms in the integrand defining u are bounded along 
C_r;-i5 u{t;x) is likewise bounded, thus implying the desired inequality to show condition (B.3). 

Condition (B.4) follows from residue calculus. In order to check it, however, we must first 
determine what initial data corresponds to step Bernoulli ASEP initial data. 

Lemma 4.15. For step Bernoulli initial data with density parameter p G (0, 1] and x G W'', 

k 



E 



.i(o)---Q.,(o) 



Lxi>0 



n 

i=i 



pr^ \pT 



k-j+i 



rj_i- 



(39) 
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with the convention that xq = 0. 

Proof. From the definition of Qx one readily sees tliat 

Qxi(O) • • • Q.,(0) = n 7?.,r('=-^^Kr('-^'+^)(^?-i-i+^ (40) 

Tliis expression involves two types of terms: r]xT^^'' and r^^"". Observe that 



E 



T 



Taking expectations of (|40p and using the above formulas we get the desired result. □ 

Thus, in order to show (B.4) we must prove that 

k 

\irau{t-x) = l,.,>onP^'"'(P^'"'^' + 1 - pT^-""^-'-'. (41) 
*^ i=i 

(Note that for p = I this simply reduces to lxi>o 11^=1 ''"^^~^)- '^^^ ^^^^^ observation is that we can 
take the limit of t — t- inside of the integral defining u{t; x). This is because the integral defining u 
is along a finite contour and the integrand is uniformly converging to its i = limiting value along 
this contour. 

When t = the exponential term in the integrand of (j37p disappears. If xi < then the integrand 
no longer has a pole at zi = — r. Since there are no other poles contained in the zi contour, Cauchy's 
theorem implies that the integral is zero, hence the condition that n(0; x) = is satisfied. 

Alternatively we must consider the case where < xi < X2 < • • • < x^. We can write n(0; x) as 

7i(0;x)=rV('=-^)/2^i(xi,...,Xfc), (42) 

where we define (for £ > 1), 

(-1)^^ f f rt ^A-ZB ^ f i + zj Y'-' 1 __e 
gi{xi,...,xk) = , / ••• / 11 11 — r ~, i^^^i- 

{27rif J J ZA - TZB V 1 + Zi/Tj T + ZiZi- T^O 

As a convention, when = we define g^ = \. 
Lemma 4.16. For i> 1 and < xi < X2 < • • • < Xk, 

5£(xi,...,xfc) = j ___g^^^(a;^,...,3;^_^). 

Proof. The lemma follows from residue calculus. Expand to infinity. Due to quadratic decay in 
Zk at infinity there is no pole. Thus, the integral is equal to —1 times the sum of the residues at 
Zfc = T~^Zj for j < k and at Zk = t^O. 

First consider the residue at Zk = T~^Zj for some j < k. That residue equals an integral with one 
fewer variable: 

(_l)fc-i r f ZA-ZB f l + z^ Y^~^ 1 dzi 



(27rO''~W J ZA-TZB \l + Zi/r t + Zi Zi 



i=i X^ + ^j/'^'^J T + Zj/r Zj/r - 
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This integrand has no pole at Zj = — r. This is because the new factor contains (1 + Zj/TY''~^ 
in the numerator and, since > Xj, this factor cancels the pole coming from the denominator 
(1 + Zj/T)^^~^. Since the contour for Zj was a small circle around — r the fact that this pole is no 
longer present implies that the entire integral is zero. This shows that the residue at Zk = zj/t for 
any j < k is zero. 

The remaining residue to consider is from = t^9. One readily checks that evaluating this residue 
leads to the desired recursion relation between ge{xi, . . . , x^) and ge+i{xi, . . . , Xk-i)- Finally, note 
that when k = 1 the recursion holds under the convention which we adopted that without any 
arguments, g£ equals 1. □ 

We may now conclude the proof of condition (B.4). Iteratively applying Lemma |4. 161 leads to 



5i(xi,...,x,)=n('^i±^') 



After some algebra one confirms that plugging this into (|42|) leads to the desired equation of (j41|) 
and hence completes the proof of condition (B.4). □ 

4.4. ASEP moment formula. We seek to compute an integral formula for the moments of Qx- 
Even if we were to solve in Proposition 14.111 (B), this would not do since x is restricted to lie in 
. Instead, the following lemma shows that we may recover the moments of Qx from the formula 
given in Corollary 14.141 for E Qx^ {t)' ' ' Qx^ (0 ■ Theorem 14.191 below gives the final formula for 

^[{Qxn 

Lemma 4.17. Recalling Qx and Qx defined in \2S^) . we have 



{QxT = jz{^^ {r-rU-l? Q^.---Q^.^ (43) 

fc=0 ^ ^ T xi<---<Xk<x 

where the empty sum (when k = 0) is defined as equal to 1. 

Proof. This lemma can be found as Proposition 3 in [19] . The derivation provided therein utilizes 
the Uq{sl2) symmetry of the spin chain representation of ASEP. We provide an elementary proof. 

Recall that Qx and Qx are functions of the occupation variables r/. If 7] is not left-finite, then 
both sides above are zero. 

In order to prove the identity, we develop generating functions for both sides and show that they 
are equal. Multiply both sides of the claimed identity by u'^/{T]T)n and sum over n > 0. The 
r-binomial theorem (see Section [7.11 with q replaced by r) implies that the generating function for 
the left-hand side of (1331) can be summed as 



For |m| small enough, this series is convergent and it represents an analytic function of u. 

Turning to the generating function for the right-hand side of (|43p . if |it| is sufficiently small, it is 
justifiable to rearrange the series in u into 

oo oo n—k 

E E (-i)Vq.,...q.,^ " 



k=0 xi<X2<---<xi^<x n>k ^ 
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The summation over n> k can be evaluated as l/(u; t)qq and factored out. Also, tlie summation 
over k and ordered sets xi < X2 < ■ ■ ■ < < x can be rewritten yielding the right-hand side of 
(jlSj) equals 

Uy<.a-uQy) 

The above manipulations are justified as long as \u\ is small enough, due to the fact that all but 
finitely many of the Qy are zero. 

The proof now reduces to showing that 

1 _ Uy<.a-uQy) 

This, however, is an immediate consequence of the definitions of and Qy. To see this, assume 
that r]y = for all y < x except when y = m , . . . , . Then Qx = and the left-hand side can be 
written as 

(l-n)---(l-ur^-i) 



(w;r)c 

,^ — n fr^y- oU n 



On the other hand, note that Qy = for all y < x except Q„. = r* ^. Thus the right-hand side can 



also be rewritten as 

{1 - u) ■ ■ ■ {1 - ut'-^] 



iu;T)oo 

hence completing the proof of the lemma. □ 

Using Corollarv 14. 141 and the symmetrization identities contained in Lemma 17.21 we can evaluate 
part of (j43p via the following result. 

Lemma 4.18. For all k > 1 

xi<-<Xk<x \ J J J \<A<B<k ^ ^ i=l * 



where the contours of integration are all C- 



-1- 



In 



Proof. The starting point for this is the formula provided in Corollary 14.141 for E Q^^ ■ • ■ Qx^ 
that formula, set = {1 + Zi)/{1 + Zi/r) and note that the contour C-r--i can be chosen to be 
a sufficiently small circle around — r so that > 1 as Zi varies in C-r--i. The summation over 
xi < • • • < Xfc < X on the left-hand side of (j4.18p can be brought into the integrand and is performed 
by using (here we rely upon > 1 for convergence) 

xi<---<Xf,<xi=l i=l ^ ^* 

After performing the summation as above, we observe that since all contours are the same, we may 
symmetrize the left-hand side. For the same reason we may symmetrize the right-hand side integrand 
in (j44p . The symmetrization is achieved by using the two combinatorial identities in Lemma 17.21 - 
identity (|67p is used to symmetrize the left-hand side, while (j68p is used to symmetrize the right-hand 
side. The two resulting symmetrized formulas are identical, thus yielding the proof. □ 



We may now prove the following moment formula. 
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^1, ^2, ^3 




Figure 5. The contour for includes 0, — r but does not include —1, t9 or 
{tzb} B>A- The dotted lines represent the images of the contours under multipli- 
cation by T. For instance, observe that the zi contour does not include the image 
under multiplication by r of Z2 or z^. 



Theorem 4.19. Fix non-negative real numbers < p < q (normalized by p + q = 1) and set all 
bond rate parameters ax = 1. Consider step Bernoulli initial data with density p £ (0, 1]. Then for 
all n > 1, 

where the integration contour for za includes 0, — r but does not include —1, t9 or {tzb}b>A (see 
Figure\^for an illustration of such contours). 



E 



Proof. The left-most equality of ()45p is just by definition. The proof of the second equality relies 
on the following lemma. For an illustration of the types of contours involved, see Figure [5j 

Lemma 4.20. Fix n > 1. Assume f{z) is a meromorphic function on C which has no poles in a 
ball around and which has /(O) = 1. Let be a small circle centered at and be another 
closed contour. Assume that there exists r > such that tC^ is not contained inside r^C^, and 
such that f has no poles inside r'^C^. Define Cf = r'C^, C} = and Q = U C} for I < i < n. 
Let 

"""TF^r/ ■■■/ n ^i^n/(-.)- 

(ZTTi)" Jc , ZA — TZb^^ Zi 
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and 

with the convention that uq = 1. Then 



-'-tit) 

k=o ^ ^ 



k(k-\) n(n-l) 

r 2 2 



Proof. In order to evaluate the integrals defining z/„ we split them into 2" integrals indexed by 
S C {1, . . . , n} which determines which integrations are along Cf (all Zi with i E S) and which are 
along Cf (all Zi with i ^ S). This shows that 



\S\=k 

where = Ij^^, 1 < z < ?i. We now claim that for any S C {1, . . . , n} with l^l = /c, 

-'Oi JO„ i<A<B<n ^ " i=l * 



(2^0" Jcji 



where we use the notation ||5|| = '^^^gi- Note that is C'^ when i G S" and C} (and hence C-^) 
when i ^ S. To prove this claim, label the elements oi S as ii < i2 < ■ ■ ■ < ik- By the fact that Zi^ 
is contained in tCj for all j > ii, we may shrink the Zi-^ contour to zero without crossing any poles 

except at = 0. The residue at that pole is t"^""*^-*. Then we may shrink the Zi^ contour to zero 
with contribution of r"^"'"*^). Repeating this up to Zj^. yields a factor of 



k 



n 

The remaining integration variables can be relabeled so as to yield the expression for Vn-k- 
By using (I46p we find that 

n 

fc=osc{i, .••,".} 

\S\=k 



fc(fc+i) .^-^ iig|| fc(fc+i) 



fc=0 Sc{l,...,n} 
|S|=fc 

n / 

_„fc+M^±i) / n 

i^n-k' 

k=0 

' n\ fc(fc-l) n(n-l) 



A;, 



k=0 



as desired. From the first line to second line is by factoring. The second line to third is by (I62p . 
The third line to fourth line is via changing k to n — k. □ 
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We return now to the proof of Theorem I4.19[ Consider the second equahty in (j45p . By virtue 
of the conditions imposed on the contours, we may apply Lemma 14.201 with f{z) = Fx-p{z) and Cj 
chosen to match the contours defined in Theorem 14.191 This shows that 

where the integration contours are C-r -i (which coincide with from Lemma l4.20p . By Lemma 
14.181 we rewrite (jITP as 

fc=0 ^ ^ T Xl< - <Xk<X 

where the empty sum (when A; = 0) is defined as equal to 1. Lemma [4.171 implies that this expression 
equals E [(Qx)"], proving the theorem. □ 



5. From nested contour integrals to Fredholm determinants for ASEP 

Using the nested contour integral formula of Theorem 14.191 for E |^r"^^^*)] under step-Bernoulli 
initial data for ASEP, we prove Mellin-Barnes and Cauchy type Fredholm determinant formulas for 
the eT--Laplace transform of r^''^*). This transform characterizes the distribution of N^it) and is 
the starting point for asymptotic analysis. The Mellin-Barnes type formula we discover is new. The 
Cauchy type formula is, after inverting the eT--Laplace transform, equivalent to Tracy and Widom's 
ASEP formula for step Bernoulli [42] initial data (see also [391 140] for step initial data where p = V). 

The route from the nested contour integral of Theorem 14.191 to the Fredholm determinants is 
similar to what was outlined in Section 13.11 (for the Mellin-Barnes type) and Section 13.21 (for the 
Cauchy type). There are, however, some differences due to the nature of the nested contours. For 
g-TASEP the integration contour for za was on a single connected contour and the set of such 
contours (as A varied) was nested so that the za contour contained {qzs} b>a- For ASEP, the 
integration contour for za is the union of two contours and the set of such contours (as A varies) is 
chosen such that the za contour does not contain {qzs} b>A- This difference in contours necessitates 
an analogous result to Proposition 13.21 (given below as Proposition 15. 2p when developing the Mellin- 
Barnes type formula, and an analogous result to Proposition 13.81 (given via the combination of 
Lemmas 14.171 and 14.181 above) when developing the Cauchy type formula. 



5.1. Mellin-Barnes type determinant. 

Definition 5.1. Fix a G C \ {0} and consider a meromorphic function f{z) which has a pole at a 
but does not have any other poles in an open neighborhood of the line segment connecting a to 0. 
For such a function and for any A; > 1, define 

where the integration contour for za contains 0, a but does not include any other poles of / or 
{TZB^ByA- For instance, when / is as in (I36p and a = — r, then the contours illustrated in Figure 
[5] suffice (for = 3) . 
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Proposition 5.2. We have that for fj.^ o,s in Definition \5.1 

1 (1-r)^ 
mi\m2\ ■ ■ ■ {2ttlY^^'^ 



Xhk 
A=l'"i2™2. 





-1 




WiT^i — Wj 





i,j=l j = l 



(49) 

where the integration contour for w j contains 0, a and no other poles of f, and it does not intersect 
its image under multiplication by any positive power of t (see Figure\B^. 

Proof. The proof is via residue calculus and follows in the same manner as Proposition 13.21 whose 
proof is found in [8] as Proposition 3.2.1. Rather than repeating that proof, we just illustrate the 
k = 2 case. 

Consider /x^ as in Definition 15.11 with contours like in Figured Initially the zi contour is chosen 
so as not to contain tz2. Because the contours must include a and 0, they must be composed of 
two disjoint closed parts. Around a the contours can be the same small circle, but around 0, the 
Z2 contour must have radius which is at least times that of the z\ contour. We wish to deform 
these contours to the single contour specified in the proposition (on which the Wj is integrated). In 
order to do this we fix the Z2 contour and deform the portion of the zi contour surrounding a so 
that it now contains the image of the circle around a on which Z2 is integrated. By the hypotheses 
on / laid out in Definition 15.11 this deformation can be performed so that the only pole which is 
encountered is when z\ crosses tz2 (we are thinking of Z2 as fixed). On account of crossing a pole, 
we find that ^2 can be expressed as 

{2Tny J J Zi- TZ2 Zl Z2 2TTi J Z2 

In the double integral it is now possible to deform (without crossing poles) the zi and Z2 contours to 
the single contour on which Wj is integrated. Similarly, the single integral contour can be deformed 
(without crossing poles) to the Wj contour as well. Observe that there are also two terms contained 
in the right-hand side of (09]) - one term is a single integral and one is a double integral. The 
single integral term matches exactly while to match the double integral we simply symmetrize the 
integrand (as can be done since zi and Z2 are on the same contour) and find those terms match as 
well. In general, the partition A indexes the clustering of residues into chains. □ 

Using the above result as well as Proposition [3?3l we find following Fredholm determinant formula 
for the Cr-Laplace transform of Qx{t) = r^''^*^ 

Theorem 5.3. Consider ASEP with < p < q (normalized by p + q = 1), all bond rate parameters 
Ox = 1, and step Bernoulli initial data with density parameter p G (0,1]. Then with notation 
9 = we have that for all ( £ C \ M+ 



E 



(CT^-W;r) 



det(/ + K^^^^^) 



where det{I + K^^^^) is the Fredholm determinant of K^^^^ : L2(Co,_r;-i,r6») — > -^^^(C'o,-t;-i,t6») 
where Cq -t--i tO 1 positively oriented contour containing 0, — r on its interior and with —1 and t6 
on its exterior. The operator is defined in terms of its integral kernel 

Kf^^{w,w') = / r(-.)r(i + .)(-c)- '-'f'^\ , ds. 
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Figure 6. The Wj contour is chosen so as to contain 0, a and no other poles of / 
(such as the one indicated with a black dot to the left of a). 



The contour D^ d is given in Definition 13.51 with d > taken to be sufficiently small and R > 
sufficiently large so that 



inf 

W,w'eCo,-T;-l,T6 



\q^w - w'\ > 0, 



and 



sup 

W,w'£Co,-T;-l,Te 



g{q''w) 



< oo. 



The function gx-p{z) is given by 

9x;p{^) = exp 
Corollary 5.4. We have that 

¥{N^{t) = m) 



(q - p)t 



Z + T 



Z + T 



(50) 



, (T™+^C;T)oodet(/ + K^^^^^)dC 



where the contour of integration encloses C, = t for < M < in and only intersects 1R+ in finitely 
many points. 

Proof. The function ^ i— t- E [l/ (^T^'=^*^;r) ] is analytic away from = t~^'^ , for integers M > 
(this is easy to see and explained in the proof of Theorem I5.3P . It follows then that even though 
the Fredholm determinant det(/ + K^^^^) is not defined for ^ S we may still integrate along 
a contour which intersects M+ only a finite number of times. We may then apply the inversion 
formula in Proposition 17.11 hence the corollary. □ 

Proof of Theorem \5.3\ . Theorem 14.191 gives a nested contour integral formula for E |^r'^^'=(*)] . Com- 
paring it with Definition 15.11 we see that /Ufc = E [r'^^"'^*^] if the contour is chosen as the one in 
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g{z)/g{qz) where 



Theorem 14.191 and if f{z) = Fx-p{z). This function can be written as f{z) 
9{z) = gx;p{z) is given in dSO]). 

We apply Proposition 15.21 yielding an expression for as in (j49p . This matches the expression 
in (jl4p up to changing g to r and sign inside the determinant. We may therefore apply Proposition 
followed by Proposition 13.61 (with q replaced by r). At a formal level, this shows that 



E 

fe>0 



Mfc— = det(/ + Kl) = det(/ + Kl) 



(51) 



the kernels and K'^ defined with respect to Fx-p{z) and gx;p{z) as above. The contour Ca in those 
propositions should be taken to be Co,-r;-i,T95 as in the hypothesis of Theorem 15.31 In applying 
Proposition 13.61 the contour Ci^2,... should be chosen to be D^ ^ with d > sufficiently small, and 
-R > sufficiently large, and the contours Ck, k > 1 should be chosen to be Dji^d^k- From the 
definition of Co,_r;-i,T9 and gx;p{z) it is easy to check that as long as |^| is sufficiently small, the 
criteria for these to be numerical equalities is satisfied. 

Since by definition Nx{t) > and r < 1, it is immediate that r^^^i^) < l. Hence, considering the 
left-hand side of ()5ip . by choosing |,^| small enough it is justifiable to interchange the summation in 
k and the expectation. By the r-Binomial theorem (see Section [7. ip we find 



E 



((l-r)^r^^W;T)^ 



det(I + J^|). 



(52) 



This equality holds for all |,^| sufficiently small. However, the right-hand side is analytic in ^ ^ 
due to Proposition 13.61 From the definition of the left-hand side 



E 



1 



((l-r)er^^W;r)^ 



E 

1=0 



P(iVx-(t) = () 

((l-r)er^;r). 



For any ^ ^ {(1 — t)~ r~ }j\/=o,i,...) within a neighborhood of ^ the infinite products are uniformly 
convergent and bounded away from zero. As a result the series is uniformly convergent in a neigh- 
borhood of any such ^ which implies that its limit is analytic. Therefore both sides of (j52p are 
analytic for ^ ^ ]R_|_ and hence by uniqueness of the analytic continuation they are equal on this set. 

The desired result for this theorem is achieved by setting ^ = (1 — t)~^C, thus completing the 
proof. □ 

5.2. Cauchy type determinant. 

Theorem 5.5. Consider ASEP with < p < q (normalized by p + q = 1), all bond rate parameters 
= 1, and step Bernoulli initial data with density parameter p G (0,1]. Then with notation 



we have that for all C £ C 



E 



(Cr^^W;r)^ 



det(/ - (K 



ASEP\ 



(C;t) 



(53) 



where det(/ — QK^^^^) is an entire function of Q and is the Fredholm determinant of K^^^P ; 
L^(C_T-;-i) — )• L^{C-r -i) defined in terms of its integral kernel 



K 



ASEP 



{w, w') 



TW — W 

with Fx-p{w) defined in i36\} . and C-r--i is a circle around — r, chosen with small enough radius so 
that —1 is not included, and nor is the image of the circle under multiplication by r (see Definition 
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Corollary 5.6. Consider ASEP with < p < q (normalized by p + q = 1), all bond rate parameters 
ax = 1, and step Bernoulli initial data with density parameter p £ (0, 1]. Then 

F{Nx{t) =m) = -r- J- / ^^^^i^ (54) 

where the integral is over a contour enclosing ^ = q~^^ /or < M < m—1. Here det{I—CK'^^~^^^^) 
is the Fredholm determinant of K'^^'^^^^ : L'^{Cr) — t- L'^(Cji) defined in terms of its integral kernel 



and €(^) = p^~^ + — 1 OL'Ti'd Cr is a circle around zero of radius R so large that the denominator 
p + q(,(,' — ^ and ^ — lp{l — r) are nonzero on and outside the contour. As a function of C; 
det(/ - ^^rH/-A5i?P) g^^^^g_ 

Proof. This follows from Theorem 15.51 (after a change of variables) and the e,-- Laplace transform 
inversion formula in Proposition 17.11 The change of variables is 

^ ^ 1+w 
1 + w/t 

Using the equivalences given in Remark 17.31 using the definition of = p/(l — p) we find that 



e-l+p(l-T) 

Similarly we find 

1 ^{r-0{r-e) q 



T(i-r) (p+gee'-e) 

The change of variables introduces an additional Jacobian factor into the new kernel which is given 

by 

-t(1-t) 

(r-e)(r-eo' 

Finally, under this change of variables, the contour C-r--i becomes as specified in the statement 
of the corollary, but with clock-wise orientation. Changing this to the standard counter-clockwise 
orientation introduces a factor of —1 into the kernel. Combining these calculations we find 



E 



1 



(CT^^(*);r)^ 



det(I - Ci?TW-ASEP) 



(C;t)c 



where asep -g -^^ statement of the corollary. 

From Proposition 17. II it follows that 



1 /■, > 'let(/-CK™-*="'^)^,_ „1 /■ det(/ - Cii:"™^ ) 



where the integral is taken over a contour enclosing ( = q~^^ for < M < m — 1, thus proving the 
corollary. □ 

Remark 5.7. For ASEP with step-Bernoulli initial data, Tracy and Widom |42] (for step initial 
data see \39\ BO] ) arrive at a very similar formula which says 

nN,M>_.,.^[ ^'^^'-f~^ f. (55) 
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where the integral is taken over a contour enclosing C = and C = Q foi' < M < m — 1. Since 
lP(A^x(i) >m)- F{N^{t) > m + 1) = ¥{N^{t) = m), it is straightforward to go from ([55]) to dH 
since 

1 1 ^1 



(C;T)mC (C;T)m+iC iC;r)m+i' 

Going in the reverse direction uses a telescoping sum and would require an a priori confirmation 
that the right-hand side of (|55p goes to zero as m goes to infinity. 



Proof of Theorem \5.5[ Let jlk be given as in ([T6]) with /(w) = Fx-p{w) defined by ([36|) and contour 
C-r--i as in Definition 14.121 Then Proposition 13. lOl and Remark 13.111 implv that 



fc>0 



det(/ + ^K) 



(56) 



where det(I + S^K) is the Fredholm determinant of 

K{w,w') = (1 - r)- 



fH 



TW — W 



We need to check that this is a numerical equality (not just formal). Because the kernel is bounded 
as w varies along C_t- it follows that K is trace-class and hence det(/ + ^i^) is an entire function of 

e 

In order to see that the left-hand side is uniformly convergent for small enough |^| we utilize the 
probabilistic interpretation for /x. By combining Lemmas 14. 171 and 14. 181 we find that 



E 



E 

fc=0 



n 



This transformation from fl^ to E l^r"^^*-*^] is upper-triangular and hence can be inverted. One 
checks that the inverse is given by 



i-lffi, = (-l)V 



^ fc(fc-l) 



j=0 



By (IMI) we find 



[-ifj,^ = (-l)V^E (1 - T^-W) • • • (1 - r^^W-^) 



E 



l)(r 



t)---{t 



NAt) 



This probabilistic interpretation of Jit implies that < 1, hence for |,^| small enough the series 
on the left-hand side of ()56p is convergent and the equality is numerical. 

By replacing ^ = — C/(l — t) and using the probabilistic interpretation for jik to justify the 
exchange of summation and expectation (assuming |C| small enough) this left-hand side series equals 



fc>0 



{-C/{l-T)f 
krl 



E 



E 



E 



(t-''-(";t)i 



fc>0 



(C;t) 



((t^-(*) 



Since we already wrote down the Fredholm determinant for this expression in (I56p . this establishes 
the claimed result of the theorem, for |(^| small enough. 
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Finally, note that 



(C;r) 



^P(iV,(t) = fe)(C;r)fc. 



oo 



fc>0 



For any C G C and any compact neighborhood Q of C, it is clear at as /c — t- oo, the product defining 
(C) T)k converges uniformly over to a finite limit. This implies that the series is likewise uniformly 
convergent in that compact neighborhood and therefore the series is analytic in a neighborhood of 
(. As C was arbitrary this implies that the left-hand side of (j53p is an entire function of We 
showed earlier that the right-hand side is entire therefore, since the two functions of C are equal for 
1^1 small enough, by the uniqueness of analytic continuations it follows that the equality holds for 
all G C, completing the proof. □ 



There are three main parameters in g-TASEP: time t, particle label n, and the repulsion strength 
q (the Oj are also present, but play a somewhat auxiliary role). On account of this, there are many 
interesting scaling limits to be explored. We will presently focus on one which involves scaling 
(7 — )• 1 and t — 7- 00, but keeping n fixed. We show that the limit of g-TASEP corresponds to 
a certain semi-discrete version of the multiplicative stochastic heat equation (and hence also the 
O'Connell-Yor semi-discrete directed polymer partition function |29j). We then introduce the limit 
of the g-TASEP free evolution equation with boundary condition and the Schrodinger equation 
with Bosonic Hamiltonian (Proposition 12.71 (B) and (C)) and show how these limits are achieved 
from the analogous statement for g-TASEP. Finally, we remark on the fact that taking a limit of 
the Mellin-Barnes type Fredholm determinant formula for the e^-Laplace transform of g-TASEP 
yields a rigorous derivation of an analogous formula for the Laplace transform of the solution to the 
semi-discrete multiplicative stochastic heat equation. 

From this semi-discrete limit it is possible to take another limit to the fully continuous (space- 
time) multiplicative stochastic heat equation [T]. The free evolution equation with boundary condi- 
tion and the Schrodinger equation with Bosonic Hamiltonian limit to the two different formulations 
of the attractive quantum delta Bose gas. 

Definition 6.1. The semi-discrete multiplicative stochastic heat equation (SHE) with initial data 
zq and drift vector a = (ai, . . . , oat) is the solution to the system of stochastic ODEs 

dz{T,n) = V z{T,n)dT + z{T,n)dBn, z{0,n) = zq^u), z{t,0) = 

where (i?i(s), . . . , i?7v(s)) are independent standard Brownian motions such that Bi has drift di, 
and we use the notation Vz(t, n) = z{t, n — 1) — z{t, n). 

There is a Feynman-Kac path integral representation for z{t, n). Let (phe a, Markov process with 
state space Z which increases by one at rate one (this is a standard Poisson jump process whose 
generator is the adjoint of V). Let £ denote the expectation with respect to this path measure on 
(j). Define the energy of (j) as the path integral through the disorder (the white noises given by dBi) 
along 4>: 



6. Appendix I: Semi-discrete directed polymers 




Also write :£',-(</') : for Ericj)) — ^. Then 



z(r, n) = f'^(")=" e^^-('^)^zo('A(0)) 



(57) 



This path integral is essentially the O'Connell-Yor semi-discrete directed polymer partition function 
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6.1. Semi-discrete limit of g-TASEP dynamics. We now show how g-TASEP rescales to the 
semi-discrete SHE. We state the result for step initial data and then provide a scaling argument 
which makes clear the correspondence for general initial data. 

Proposition 6.2. Consider q-TASEP started from step initial data and scaled according to 

q = e-% ai = e-'^\ t = e~^T, x„(t) = e'V - (n - l)e~^ log - e"^Ff (r). (58) 

Let z^{T,n) = exp (— ^ + F^i^))- Then for any N > 1, T > 0, as e — )• 0, the law of the stochastic 
process {ze{T,n) : r € [0,T], 1 < n < N} converges in the topology of measures on C([0, r],]R^) to a 
limit given by the law o/{z(r, n) : r G [0,T], 1 <n< N} where z{T,n) solves the semi-discrete SHE 
with drift vector a = (ai, . . . ,0,^) CLnd initial data Z{){n) = 5n=i- 

This result is a corollary of [8] Theorem 4.1.26 which deals with a larger 2-dimensional extension 
of (7-TASEP and its limit. That result is not entirely elementary as it relies upon the convergence of 
g'-Whittaker processes to Whittaker processes [8] as well as the relationship of Whittaker processes 
to the semi-discrete directed polymer [28]. We will presently provide a purely probabilistic sketch 
of why this result is true, without making any attempt to fill in the details of rigorous justifications. 

It is easy to check the initial data. Observe that via the scalings, Ze(0, n) = £'""^6*^". Hence if 
n > 1 the limit is 0, whereas for n = 1 the limit is 1. This shows that as e — )• 0, ^^^(O, n) — t- 6n=i- To 
achieve a general initial data zo, one should scale Xn{0) so that g"-ie-«^i(o) — ). zo{n). 

To see how the dynamics behave under scaling it is easiest to work in terms of F^{t). Observe 
that 

dF^ir) = F:{r)-F:{T-dT) 

= (e~V- (n- l)loge~^ -ex„(e~V)) - (e~^(r-(ir) - (n - 1) log e"^ - ex„(e"V - e"^ 
= e~^dT — e (2;„(e~^r) — x„(e~^r — e~^(ir)) . 

The jump rate for q'-TASEP, in the rescaled variables, is given by 

This means that in an increment of time e~'^dT, we should see that 

e (x„(e-V) - - e~'dT)) = - (a„ + e^""'^-^"^) dr + dW^ + o(l), 

where the Wn are independent Brownian motions which arise from the approximation of a Poisson 
process by a Brownian motion. Setting Bn = dn — Wn (a Brownian motion with drift a.„ now) we 
find that, 

dF,"(T) = e^""M-^"M + dBn + 0(1). 

By Ito's lemma, 

dexp (FJ^(t)) = Q exp (Ff (r)) + exp (F,"" \r))) dr + exp (F,"(r)) dB^ + o(l), 

and hence rewriting this in terms of z^ (r, n) we have 

dZf:{T, n) = Vze{T, n)dT + z^{t, n)dBn + o(l). 
As e — )■ this equation limits to that for z(t, n) as desired. 
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6.2. Semi-discrete limit of g-TASEP duality. By utilizing the path integral formulation of 
z{t, n) given in (j57p let us compute expressions for joint moments of z{t, n) for fixed r but different 
values of n. For simplicity, we assume below that all = 0, though the general case is no more 
difficult. This procedure is sometimes called the replica approach (not to be confused with the 
replica trick - see Section 16. 4p as it involves replication of the path measure. 
Observe that 



E 



.i=l 



E 



JJ^</'i('r)="i 



.i=l 



where the (p'^s are independent copies of the Poisson jump process (p. Interchanging the disorder 
and path expectations we are left to evaluate the (now inner) expectation 



E 



.i=l 



■■Er{(p^)■. 



exp 




This leads to the final formula 

k 



E 



Y[z{t, Ui) 



.1=1 



g'l>l(T)=ni _ _ _ g(j>k{r)=nk 



exp 



<i>i{s)=<f>j{s) 



n^o(<^i(o)) 



(59) 



This identity should be thought of as a duality between the semi-discrete SHE and a system of 
Poisson jump processes energetically rewarded via the sum of their local times. The proof of the 
above identity follows from the simple fact that for X distributed as a centered normal random 
variable with variance o"^, 



E 



This implies that it is the Gaussian nature of the noise and not the underlying generator V which 
is behind this identity. Therefore, if V is replaced in Definition 16 . 1 1 by an arbitrary generator L, the 
same identity holds if (j) is defined via the adjoint generator of L. For more on these generalities, 
see Section 6 of [8]. Note that for the continuum SHE, there exist other types of noise for which 
dualities have been shown (see for instance |18j). 

Just as for (7-TASEP, (j59p implies that the joint moments of z satisfy systems of ODEs (recall 
Proposition 12. 7p . The (A) system follows from (I59p directly. We now record the limiting versions 
of Proposition [221(B) and (C). 



Proposition 6.3. Let z{T;n) be as above with initial data ZQ{n) supported on Z>o. 

) Free evolution 
(1) For all n G 



(B) Free evolution equation with boundary condition: If -u : IR+ X (Z>o)^' M solves: 



and T G 



d 



■^u{t; n) = ^ Viu{T; n); 



i=l 



(2) For all n G 



such that for some i G {1, . . . ,k — 1}, rij = "/ij+i, 



(3) For all n G (Z>o)^ such that Uk = 0, u{t; fi 



(4) For all n G W^^, u(0; 



n 



(Vi - Vi+i 
n«=i^o(rai); 



1)u{t; n) = 0; 

for all T G 
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Then for all n G W^^r., E 



nLi2:(T,ni) =u{T;n). 



(C) Schrddinger equation with Bosonic Hamiltonian: If v : x (Z>o)'^ solves: 



[1) For all n E (Z>o)'' and r G M+, 



—v{T;n) = Hv(T;n), H 
dr 



i=l 



(2) For all permutations of indices a G Sk, v{T;aii) = v{T]n); 

(3) For all n G (Z>o)^ siic/i t/iat = 0, v{t; n) = /or a// r G M+; 

(4) For all n G VF^q. *(0;n) = [iti ^o(ni); 



T/ien /or all n G W^^n, E 



Yti=iz{T,ni) =v{T;n). 



These systems of ODEs can be proved from (|59|) directly. Instead we sketch how they arise as 
hmits of the analogous ODEs for g-TASEP. 

Let us first consider (B). Recah that u{t;n) = E OLi 9^^"'''*^'^"'' • Thus, defining 



u, 



(r,n) = ]Je^"'"e("'~i)u(e^V, 



n) 



i=l 



we expect (from Section 16. 1|) that 



_3fcr 

lim e 2 u^(t, n) = E 



n 

.i=l 



Call this limit u{T,n). We now check that u indeed satisfies conditions (B.1-4) above. The fact that 
it satisfies (B.3) and (B.4) is clear. Note that 



(60) 



i=l 



Using this, it follows by rescaling (B.l) of Proposition [2771 that 
d 



— Ue(r, n) = ke ^Ueir, n) + (e ^ - |) ^ (e'«e('r, n- ) - n,(T, n)) + 0(e). 



1=1 



The factor e ^ ^ comes from the expansion of e {1 — q). The above can be rewritten as 

d 



i=l 



which in turn implies that 



3kT 



— e 2 n,(r,n) = 2^Vie 2 n,(r, n) + 0(e) 



1=1 



This shows that in the e — )• limit, u satisfies (B.l) above. 

Using (j60p and the expansion g = 1 — e + O(e^), it follows from (B.2) of Proposition 12.71 that 
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__3fcT 

Multiplying by e 2 has no effect on this equality, and so in the limit e — )■ 0, we find that u satisfies 
(B.2). 

We now consider (C). Define and v analogously to and u above. The fact that v satisfies 
(C.2), (C.3) and (C.4) is clear. Using (j60p and second order expansions of (1 — q) and (1 — q~^) we 
find that 

d ^ 

— i;,(r, n)=^ (i;,(r, n~) + it;,(r, n)) - ^ 5n,=n,Ve{T, n) + 0(e). 

_3fcr 

Multiplying by e 2 and taking e — )■ leads to (C.l) as desired. 

For zo(ra) = J][*L]^ 5ni=i initial data, it is possible to explicitly solve (B) and (C) in Proposition 
16.31 via nested contour integral formulas which arise as scaling limits d?]). In fact, if we change the 
boundary condition in (B.2) to (Vj — Vj+i — c) for any c G M (or analogously put this c factor in 
(C.l) in front of the sum over i < j) essentially the same integral formulas work and we find that 
(B) is solved by 

= /o \k I I n — n -"T^^i (61) 

^ ' \<A<B<k ^ ^ j=i J 

where the integration contour for wa contains and {wb + c}b>a- These systems of ODEs are 
semi-discrete versions of the delta Bose gas, and c plays the role of the coupling constant. This 
remarkable symmetry between attractive (c > 0) and repulsive (c < 0) systems is discussed more in 
Section 6 of [8]. 



6.3. Semi-discrete limit of (7-TASEP Fredholm determinant. Proposition 16.21 implies that 
as (/ — 7- 1, under proper scaling g-TASEP converges to the solution of the semi-discrete SHE. From 
this weak convergence result it follows that the e^-Laplace transform of particle location for q- 
TASEP converges to the Laplace transform of the limiting SHE. This Laplace transform completely 
characterizes the one-point distribution of the solution z(t, n). The g-TASEP Mellin-Barnes type 
Fredholm determinant formula has a nice scaling limit and thus yields (we will state it for a zero 
drift vector): 

Theorem 6.4. For r G M+, and n > 1, the solution of the SHE with delta initial data and drift 
vector a = (0, . . . , 0) is characterized by (for Reu > 0): 



E 



3t 

g-?ie"2"2(T,n) 



det(/ + Ku) 



where det{I + K^) is the Fredholm determinant of Ku '■ L^(Co) — )• L?'{Cq) for Cq a positively oriented 
contour containing zero and such that for all v,v' £ Cq, we have \v — v'\ < 1/2. The operator Ku is 
defined in terms of its integral kernel 



1 /■/,oo+l/2 

Kuiv,v') = — dsri-s)Til + s)-. 

J-(,oo+l/2 



-.00+1/2 V(s + 7;)" v + s-v' ' 



Proof. This is proved in [8], Theorem 5.2.10. An alternative choice of contours is developed in [9] 
Theorem 1.16. The formula follows from rigorous asymptotic analysis of Theorem 13. 121 □ 
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6.4. The replica trick. It is enticing to think that one might be able to compute the Laplace 
transform formula in Theorem 16.41 directly from the explicit formula for K[z{t, n)''] (such as the one 
given by combining (j6ip with Proposition 16.31 (B)). If X is a suitably nice non-negative random 
variable (for instance, if X were bounded), then for u with Re(n) > 0, 

k=0 

This identity only makes sense if one can rigorously justify interchanging the summation. Yet 
worse, if the moments of X grow too rapidly, the right-hand side might not even be convergent 
for any value of u even though the left-hand side would be necessarily finite. This is exactly the 

3r 2 

case when X = e 2 z(t, n). From ()6ip one can estimate that for this choice of X, E[X ] ~ e'^*^ 
where Cfc > c > for all k. This means that, from a mathematical perspective, one can not use this 
approach to compute the Laplace transform. 

One variation of the so-called replica trick discussed in physics literature is an attempt to sum 
this divergent series in such a way as to guess the Laplace transform. (In fact, the most typical 
version of the replica trick asks for less than the Laplace transform, rather just for E[log z(t, n)], 
and tries to access it from analytically continuing formulas for integer moments to k = 0.) 

This replica trick procedure has been implemented for the continuum SHE (a scaling limit of the 
semi-discrete SHE) in which the ODEs in Proposition 16.31 (B) and (C) become two equivalent forms 
of the attractive quantum delta Bose gas. The moments of the solutions of the continuum SHE 
grow even faster, like e'^*='^^ for > c > 0. By diagonalizing the Bosonic Hamiltonian (the limit of 
(C)) via the Bethe Ansatz, \15\ both made initial attempts at computing the Laplace transform 
via the replica trick. These initial attempts yielded a wrong answer. However, very soon afterward, 
the formula of [Il[33] was posted (with a rigorous proof given in [1]) and [El [11] showed that their 
approach was able to recover the correct Laplace transform formula. 

7. Appendix II: Combinatorics 

7.1. Useful (/-deformations. We record some (/-deformations of classical functions and transforms. 
Section 10 of p^j is a good reference for many of these definitions and statements. We assume 
throughout that \q\ < 1. The classical functions are recovered in the q ^ 1 limit. 

The q-Pochhammer symbol is written as (a; q)n and defined via the product (infinite convergent 
product for n = oo) 

(a; q)n = (1 - a)(l - aq){l - aq^) • • • (1 - ag""^), (a; q)oo = (1 - a)(l - aq){l - aq^) ■■■ . 

There are two different q-exponential functions which were introduced by Hahn [TTj in 1949. The 
first (which we will use) is denoted eq{x) and defined as 

while the second is defined as 

Eg{x) = ( - (1 - q)x;q)^. 

Both eq(x) and Eq{x) converge to as gr — )■ 1, cf. ([63]) below. In fact, eq{x) converges uniformly 
to on X S [— oo, a] for any a S M. 

The q-factorial is written as either [n]q\ or just nq\ and is defined as 

iq;q)n (1 - (/)(1 - g^) • • • (1 - (/") 

71 ' = = 

{1-qr {l-q){l-q)---{l-q) ' 
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The q-binomial coefficients are defined in terms of g^-factorials as 



n 



iQ;Q)r 



kj g kq\{n-k)q\ {q;q)kiq;q) 



n—k 



We also have |22 



= E 

1 5e{l,...,n} 
\s\=k 



||5||_M^ 



(62) 



where 

\\S\\=^i- 
ies 

The q-binomial theorem ([2] Theorem 10.2.1) says that for all |x| < 1 and \q\ < 1, 

_ («^;9)oo 



k=0 



{q;q)k {x;q)c 



Two corollaries of this theorem ([2] Corollary 10.2.2a/b) are that under the same hypothesis on x 
and q, 

oo fc ( fc — 1 ) 



OO i. 



■{-xy 



k=0 



k=0 



kq\ 



For any x and g we also have ([2] Corollary 10.2.2.c) 

Define the following transform of a function / G ^^(Z>o): 



(63) 



(64) 



n>0 



f{n) 

(C9";9)c 



(65) 



where C £ C. 



We call this the e„-Laplace transform of q-^ since if X is a random variable taking values in 



^>o 



and f{n) = P(X = n) 



no = E 



Cq 



X 



l-q 



An inversion formula is given as Proposition 3.1.1 of [8j and can also be found in [3]. 

Proposition 7.1. One may recover the function f E ^^(Z>o) from its transform /"^(C) with ^ S 
C \ {9~'^}fe>o I'iC' the inversion formula 

/(m) = -(z'"^/ {q'^^'CqUfiOdC, (66) 

where Cm is any positively oriented contour which encircles C = Q~^^ for < Ad < m. 
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7.2. Symmetrization identities. We state and prove the following two useful symmetrization 
identities. 



Lemma 7.2. For all k > 1 



- (t; r)fcT 2 zi---Zk det 



E n 

^S,l<A<B<k MA)-rz.iB) 



Zi TZ^ 



k 

«J=1 



(67) 



Setting 



l+z. 



1+^ 



4- we aZso have 



T. n 



n 



fc(fc-i) 



(-l)V 2— det 



Zi TZ-j 



i,j=l 1=1 



Remark 7.3. Before proving these identities note that for = y-^z )t 



Zi Zj 



q 



- 6 



(zi + l){p + qzi) 

T + Zi 



p + qCi^j - 
pC^ + q^i - 1> 

r- 1 



(69) 



1 - ^^/r 



Proof of Lemma \7.^ The first identity ([67|) is [25], 111(1.4). The second identity is equivalent to 
the identity (1.7) in [39]. In order to see this equivalence expand the Cauchy determinant as 



det 



fc(fc-i) 
r 2 



=1 ^i---^;fc(i-T)'' 



n 



l<i^j<k 



Zi TZi 



Multiply both sides of the claimed identity by the factor rii<j-^j<fc -j-^^, reducing the identity to 



E n 



o-GSfe l<A<B<k 



1 Ziil 



Za{B) - Z„(A) fj^ ' ' ' "1 fj[ Zi{l - t) 

Noting that "l^'^^V = (^j — 1)"^ and using the relation (1691) it remains to prove that 

P + q + ^a(B)^a{A) - ^a{A) -TT 1 



fc(fc-l) 

2 



n 

l<A<B<k 



n 



Using the anti-symmetry of the Vandermonde determinant we rewrite the above as 

1 _ k(k-i) llA<Bi^B - U) 



T=n(f.-i)-'- 



J] Sgn(o-) Yl iP + q + LiB)^a[A)-CaiA))Yl-r- t _1 

ae5fe \<A<B<k i=\ ^<^(1) ■ ■ ■ ^^(*) 

The above identity is (1.7) in [39], and the proof is complete. 



□ 
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7.3. Defining a Fredholm determinant. Fix a Hilbert space L^(X, /^t) where X is a measure 
space and /i is a measure on X. When X = T , a. simple smooth contour in C, we write L'^iT) 
where (j. is understood to be the path measure along T divided by 27rt. When X is the product of 
a discrete set D and a contour F, /i is understood to be the product of the counting measure on D 
and the path measure along F divided by 2-Ki. 

Let K be an integral operator acting on /(•) G L'^{X,fi) by {Kf){x) = J-^ K{x,y)f{y)d/j,{y). 
K{x, y) is called the kernel of K. A formal Fredholm determinant expansion of / + i^T is a formal 
series written as 

oo ^ „ „ n 

det{I + K) = l + Y.- / det[K{x,,x,)]l^^^l[dfi{xi). 

n=l ^' •'^ ''^ ' 1=1 

If the above series is absolutely convergent then we call this a numerical Fredholm determinant 
expansion as it actually takes a numerical value, li K is a trace-class operator (see [23] or [7]) then 
the expansion is always absolutely convergent, though it is possible to have operators which are not 
trace-class, for which convergence still holds. 



8. Appendix III: Uniqueness of systems of ODEs 

We prove the uniqueness result of Proposition 14.91 by a probabilistic approach. It is possible to 
extend this proof to a more general class of generators, but we do not pursue this here. 



Proof of Proposition \4-y[ Let us first demonstrate the existence of one solution to the system of 



ODEs given in Definition 14.81 Let denote the proposed solution, equation ()3ip . in the statement 
of the proposition. The definition of the generator implies that satisfies condition 1 of Definition 



To prove that satisfies conditions 2 and 3 requires an estimate. In time t, the number of jumps 
in ASEP is bounded by a Poisson random variable with parameter given by constant time t. This 
means that for some constant c' > 



P(p(-t) - x||i =n)<e 

Observe now that 

h^{t;x) - p-*'^(f(0) = x)hoix) < 



-c't 



(c't)'' 



(70) 



< 



E E 

n>l x':\\x—x'\\=n 



p-*'-MO) = f')|/io(x)l 



-Ce- 



n>l x' 



X' : \\x—x' \\=n 



n>l 

< e^ll"1li(e'="'*-l) 

The first inequality follows from the definition of as an expectation, along with the triangle 
inequality. For the second inequality we can use the bounds jSO]) and jTOD to replace p-*;^(f (0) = 
x')|/io(x')| by e~'^'*^^^^:^Ce~'^ll^'lli. For the third inequality we observe that \\x'\\i > max(0,n— 
Plugging this bound into e"'^"^ we find that the summand is now independent of x' and the 
summation over x' can be replaced by a rough combinatorial bound of (c")"" for the number of such 
x' (c" is some sufficiently large constant) - this provides the fourth inequality. The final equality 



FROM DUALITY TO DETERMINANTS FOR Q-TASEP AND ASEP 



48 



comes from factoring out e'^H^H^ from the summation and then bounding the remaining summation 
in ?i > 1 by the Taylor series for the exponential. 

The conclusion of the above line of inequalities is that for some c'" > 0, 



h\t;x) - p-*;^(x(0) = x)ho{x) 



!)• 



Observe that using the triangle inequality and the exponential bound on /io(x), the above inequality 
implies that h} satisfies condition 2. Similarly, as t — ^ 0, P^*'^(x(0) = x) — t- 1 and 6*^11^11^ (e'^'"* — 
1) — 7- we obtain the pointwise convergence (condition 3): 

h^{t;x) ho{x), t — ;> 0. 

The argument to prove uniqueness is very similar to the argument used to prove condition 3. 
Assume now that in addition to h^, there existed another solution to the true evolution equation 
which we will denote by h?. The idea is to prove that g := h} — must be identically 0. The 
solution g has zero initial data. 

To prove that (7 = it suffices to show that for any T > and any x G , g{t',x) = for all 
t E [0,T]. We readily see that for any 5 G (0,T], 



g{t;x) 



n{T) 

E E 

".=0 x':\\x-x' 



< 



< 



p-*'^(f(-5) = x')g{5;x') 

=n 



E ^ E 

n>n{T) x':\\x—x'\\=n 
n>n{T) ^' 

where n(T) is a positive integer which depends on T and will be specified soon. The above inequal- 
ities follow for similar reasons as in the proof of condition 3 for h^. Now observe that by choosing 
n(T) sufficiently large, the summation in the last line above can be made arbitrarily small. This is 
due to the fact that 1/n! decays super-exponentially. This means that for any e > and any T > 0, 
there exists n(T) such that 

n{T) 

g{t;x)-Y^ Yl 1P-*'"(^M) =x')5(5;x') <e. 

n=0 x':\\x—x'\\=n 

Since the set of x' such that — x'|| = n with n G {0, 1, . . . , n{T)} is a finite set, condition 3 implies 
that as 5 —7- 0, each g{5; x') — )■ as well. Choosing 6 sufficiently small, this implies that 

\9{t;x)\ < 2e, 

and since e was arbitrary this implies in fact that g{t;x) = for all t G [0, T]. This completes the 
proof of uniqueness. □ 

9. Appendix IV: GUE Tracy- Widom asymptotics for ASEP 

We provide a critical point analysis for the long-time asymptotics of our Mellin-Barnes type 
Fredholm determinant formula for the e^-Laplace transform of t^'^^^^ We assume that r < 1 is 
fixed and straightforwardly arrive at the GUE Tracy- Widom limit theorem recorded in ([T]) and 
proved first (via an analysis of the Cauchy type formula) by Tracy and Widom [41]. In order to 
make this analysis rigorous one would need to control the tails of the integrand defining the kernel. 
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Another scaling limit of interest is the weakly asymmetric limit in which r goes to 1 simultaneously 
with t going to infinity. Under the correct scaling (as in [1]) our formula should lead to the Laplace 
transform of the Hopf-Cole solution to the KPZ equation with narrow wedge initial data. We do 
not pursue these directions presently, but rather remark that it appears that the Mellin-Barnes type 
formula is very well suited for such a rigorous asymptotic analysis. 

For simplicity, let us consider ASEP with step initial data and fix x = 0. We seek to study the 
large t behavior of Nq^I) via its e,-- Laplace transform. Let us recall the formula we have proved in 
Theorem 15.31 



E 



e,(Cr^-W) =det{I-K^) (71) 

where det{I—K(^) is the Fredholm determinant of the operator : L^(Co,-r;-i,r6») — ^ L'^{Co-T;~i,Te) 
defined in terms of its integral kernel 

expt^ ds 



Kc{w,w') = ^ [ T{-s)r{l + s)[-{l-T)CY 



expt^qiT^ 

As noted in Section [7.H er{z) converges uniformly for z G [— oo,0]. This means that if z — )• — oo 
then er{z) — )• and if z — >• then er{z) — )• 1. On account of this, if we set 

then it follows (cf. [8] Lemma 4.1.39) that 



hm E 

t—>-oo 



_^^Afo{t/7)) 



tV3 



> -r 



where we have set j = q — p- 

Theorem 15.31 (see equation (|7ip above) reduces this to a problem in asymptotic analysis. We 
proceed now without careful estimates and only discuss contours briefly. There are a few estimates 
which would be necessary to turn this into a rigorous proof. Making the change of variables in (|7ip 
z = T^w and using the fact that r(— s)r(l + s) = vr/ sin(— vrs) we arrive at the following 

t 



lim P ^ > -r^ = lim det(/ - K'A 



where the kernel is now given by 

1 1 r vr 
27rtlogr J sin(7r(log^ w — log^ z)) 



K'^{w,w') = ^t:^ I —TZTTZT^. — , (i _ ^)log, ^-log, ^ 



exp(t[G{z) -G{w)]+t^^^logTr[log^z-log^w]^ ^ 



with 



4 T + z 

The critical point of G{z) is readily calculated by solving 

G'(^) = -^ + ^^-2 =0. 

Az (r + z)^ 

This yields Zc = t as the critical point. Actually, it is a double root of the above equation and 
accordingly one sees that G"{zc) = 0. The fact that the third derivative is non-zero (and the second 
derivative is) indicates t^^^ scaling. Up to high order terms in {z — r) and {w — t), we have 



FROM DUALITY TO DETERMINANTS FOR Q-TASEP AND ASEP 



50 



The w contour can freely be deformed to go through the critical point r and to depart it at angles 
ivr/S (oriented with increasing imaginary part). Likewise the z contour can go through r — t^/^ 
and depart at angles ib27r/3 (oriented with decreasing imaginary part - as is a consequence of the 
change of variables). The w contour needs to cross the negative real axis between —1 and — r. The 
z contour must stay inside the w contour. There is one nuance with the z curve which is that it 
keeps wrapping around in a circle (since the imaginary part of s went from — oo to oo). However, 
with a suitable a priori bound one should be able to show that this s contour can be made finite at 
a cost going to as t goes to infinity. 

The above considerations suggest scaling into a window of size t^^'^ around the critical point r. 
Consider the change of variables z — t = t~^/^Tz and likewise for w and w' . This leads to: 

_J: ^ ^ ^1/3 1 



log T sin(7r log^ w — log^ z) w — z 

-,3 „..3 

t[G{z)-G{w)] ^ + — 

^ ^ ' ^ 48 48 

t^/^ log rr (log^ z — log Tw;) w r{z — w) 

(1 _ 7-)l°gr^-logr»" ~ I 

1 -1/3 



t 

z — w' z — w' 

- ^ t-'/^dz. 

z 

Additionally there is an extra factor of r^^t^^^ which is factored into the kernel, due to the Jacobian 
of the w and w' change of variables. Combining all of these factors, we get that the kernel has rescaled 
to 



^ ^ -i + €+rrz-n,) 1 dz 



2TTi J W — Z Z — w' 

where the w contour is given by two infinite rays departing 1 at angles ibvr / 3 (oriented with increasing 
imaginary part) and the z contour is given by two infinite rays departing at angles ib27r/3 (oriented 
with decreasing imaginary part). Recalling z = 2'^/^z and likewise for w and w' yields 



dz 



2-KL J W — Z Z — w' 

The Fredholm determinant with this kernel is readily shown to be equivalent to the Fredholm deter- 
minant of the Airy kernel (see for instance Lemma 8.6 of [9]), and hence its Fredholm determinant 
is equal to FGUE(2^/^r). 



This implies that 



as we desired to show. 



hm P 
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